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Experimental Elucidation of the Life Cycle of Rhinebothrium 
urobatidium (Cestoda: Rhinebothriidea) from the Round Stingray 
(Urobatis halleri: Myliobatiformes) to First and 
Second Intermediate Hosts 


Ralph G. Appy,!* Shana K. Goffredi,* Bruno Pernet,* and Carin Latino! 


' Cabrillo Marine Aquarium, 3720 Stephen M. White Drive, San Pedro, CA 90731 
* Occidental College, Dept of Biology, 1600 Campus Rd, Los Angeles, 90041 
> Department of Biological Sciences, California State University, Long Beach, 1250 
Bellflower Blvd, Long Beach, CA 90840 


Abstract.—The life cycle of the cestode Rhinebothrium urobatidium, whose final host 
is the round stingray, Urobatis halleri, includes a copepod as the first intermediate ~ 
host and small benthic fishes as second intermediate hosts. Hexacanth embryos within 
nonoperculate, untanned eggs collected from round stingray developed in the tidepool 
copepod Tigriopus californicus into caudate procercoids with an apical organ/sucker 
and cercomer, and rarely, bothridia. The procercoids, which developed in the body 
cavity within a membrane, adhered to the copepod intestine, were infective to arrow 
gobies, Clevelandia ios, within 15 d of exposure at 21°C. When infected copepods were 
fed to arrow gobies, procercoids developed into nonlacunate plerocercoids each bear- 
ing individually retractable bothridia and an apical sucker. Within 10 d of exposure 
to the second intermediate host, the larvae had migrated up the bile/cystic duct into 
the gall bladder, where they developed bothridia similar to those of adult worms. Be- 
tween 30 to 51 d post-infection in the goby, plerocercoids approached the size of larvae 
found in natural infections, and the scolex became morphologically similar to that of 
adult worms from round stingrays. Only presumptive filiform microtriches (filitriches) 
were present on procercoids while both filitrichs and spiniform microtriches (spini- 
trichs) were present on plerocercoids. Identification of plerocercoids from experimen- 
tal infections as those of R. urobatidium was confirmed through morphology of the 
scolex and using cytochrome c oxidase I sequences. The experimental transmission of 
R. urobatidium to first and second intermediate hosts provides improved understand- 
ing of the transmission and ontogeny of shark tapeworms. The biological character- 
istics of U. halleri, with its diverse parasite fauna, provide significant opportunities to 
examine the biology of an array of elasmobranch tapeworm taxa. 


The Rhinebothriidea comprises approximately 24 genera and 143 species of eucestode 
tapeworms found as adults in the spiral intestine of batoid elasmobranchs (Caira and 
Jensen 2014; Ruhnke et al. 2017; Coleman et al. 2019a; Coleman et al. 2019b). While there 
are numerous reports of larval stages (metacestodes) of members of the Rhinebothriidea 
from various invertebrate and fish intermediate hosts (Chambers et al. 2000; Caira and 
Reyda 2005; Jensen and Bullard 2010), there are no publications that describe the exper- 
imental transmission of these parasites from egg to adult worm (Caira and Jensen 2014). 


* Corresponding author: r.appy@cox.net 
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This lack of information on the Rhinebothriidea and marine cestode life cycles in general 
is surprising considering the abundance and diversity of these worms. The sparsity of infor- 
mation on marine cestode life cycles is an impediment to a more complete understanding 
of tapeworm ontogeny and evolution (Beveridge 2001) and to a better understanding of 
the important role parasites play in the food chain of many ecosystems (e.g., see Lafferty 
et al. 2006). 

During studies on fish and macroinvertebrate parasites of Anaheim Bay in Southern Cal- 
ifornia, plerocercoids with a scolex resembling that of R. urobatidium (Young, 1955) (see 
Appy and Dailey 1977) were found in the gall bladder and bile duct of small benthic fish. 
Larvae were particularly abundant in arrow gobies, Clevelandia ios, (Jordon and Gilbert, 
1882) and shadow gobies, Quietula y-cauda (Jenkins and Evermann, 1889) collected in four 
man-made tidal ponds. Rhinebothrium urobatidium is also one of the most abundant para- 
sites of the round stingray, Urobatis halleri (Cooper, 1863), which is seasonally present in 
large numbers in the same tidal ponds (Jirik and Lowe 2010) where we observed plerocer- 
coids in small benthic fishes. The purpose of this study is to describe the morphogenesis of 
the tapeworm in its intermediate hosts, and to confirm the identity of these plerocercoids as 
R. urobatidium with morphological features and DNA sequence data. 


Materials and Methods 


Round stingrays were collected at a number of collecting events from 2016 to 2018 in 
Anaheim Bay (hook and line, beach seine), Seal Beach (hook and line), San Diego Bay 
(hook and line), Long Beach Harbor (otter trawl), and Two Harbors, Catalina Island (by 
spear), California. A survey of fish parasites was conducted in Anaheim Bay in 2013 using 
an otter trawl, beach seine and minnow traps. In addition, small benthic fish examined 
for infection with R. uwrobatidium were collected in Anaheim Bay tidal ponds by beach 
seine in the fall of 2018. Arrow and cheekspot gobies, //ypnus gilberti (Eigenmann and 
Eigenmann, 1889) (Gobiidae) used for experimental infections were collected in Cabrillo 
Salt Marsh, Los Angeles Harbor. For morphological and sequence data comparison with 
larvae obtained experimentally, adult R. urobatidium were collected from round stingrays 
at multiple localities in Southern California (see above) and plerocercoids were from gobies 
collected in Anaheim Bay. 

Gravid proglottids of R. urobatidium taken from round stingrays were placed in seawater 
in small (25 ml) plastic cups, and, following expulsion of eggs, were removed from the cups. 
Eggs were rinsed with several changes of seawater. Eggs were examined with a compound 
microscope to confirm their morphology as that of R. urobatidium. Eggs were maintained 
in seawater at 21°C with daily seawater changes until use in experimental infections, which 
usually occurred within 3 d. 

Experimental infections with R. urobatidium were carried out using the tidepool cope- 
pod Tigriopus californicus (Baker, 1912) (Harpacticoidea) obtained from a colony main- 
tained at the Cabrillo Marine Aquarium as a first intermediate host, and wild caught ar- 
row and cheekspot gobies as second intermediate hosts. None of 23 gobies examined at 
the time of collection in Cabrillo Salt Marsh were infected with R. urobatidium. At a sub- 
sequent collection, | of 31 gobies was infected with a single R. urobatidium. In addition, 
infections were attempted using ostracods [tentatively identified as Cyprides beaconensis 
(Leroy, 1943)] collected in tidal pools in San Diego Bay as a first intermediate host and 
commercially available mosquitofish [Gambusia affinis (Baird and Girard, 1854)] as a sec- 
ond intermediate host. 
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Eggs of R. urobatidium were pipetted into 25 ml plastic cups containing copepods or 
ostracods. After 8 hours, copepods and ostracods were placed in a loosely covered 125 ml 
plastic cup in an incubator at 21°C. Copepods were fed Nanno 3600 (Instant Algae) daily 
and ostracods were fed finely ground marine fish pellets once a week. Evaporation from 
the cups was compensated by adding distilled water. Five to ten copepods were examined 
for larvae of R. urobatidium at each sampling event over a 20 d period. Ostracods were 
examined 15 d post exposure (dpe) to eggs of R. urobatidium. 

At 15 dpe two to five copepods experimentally infected with R. urobatidium were fed to 
both gobies and mosquitofish, that had been placed individually in 125 ml plastic cups. 
After eight hours, all copepods had been consumed and gobies and mosquitofish were 
removed from the individual cups and each species placed in a liter glass bowl at 21°C. 
Gobies and mosquitofish were fed ground dry marine aquarium fish pellets and frozen 
copepods. Gobies were examined at irregular intervals or at times when a fish was found 
moribund such as occurred at 51 dpe. Both arrow and cheekspot gobies were used in ex- 
perimental infections although measurement data presented herein came only from plero- 
cercoids removed from arrow gobies. All fish examined in the lab were euthanized using 
~250 mg/1 of Tricaine (MS222). 

Larvae isolated from copepods and dissected from gobies were observed live under a 
petroleum jelly-rimmed coverslip and then heat killed using the flame of a match prior to 
measuring. Measurements of plerocercoids were made on worms heat-killed prior to plac- 
ing the coverslip to avoid effects of compression by the coverslip and attachment/flattening 
of bothridia adhering to the slide/coverslip. Measurements, drawings, and photographs 
were made using a Wild M12 compound microscope fitted with an optical graticule, draw- 
ing tube, and a Canon EOS 60e digital camera. Hand drawings were scanned and inked 
using Adobe Illustrator. The contrast and brightness of photographs and micrographs were 
sometimes enhanced using Adobe Photoshop. Unless stated otherwise, measurements are 
given in micrometers (m). 

Selected specimens used for scanning electron microscopy (SEM) included: adult worms 
from round stingray, plerocercoids from naturally and experimentally infected arrow goby, 
and procercoids from experimentally infected T. californicus. Specimens were placed di- 
rectly into fixative or were heat killed with steaming hot water, fixed in 1.25% glutaralde- 
hyde, post-fixed in 1% osmium tetroxide, dehydrated in an ethanol series, critical point 
dried using C02, and mounted on stubs. Specimens were coated with gold/palladium us- 
ing an Emitech K550x sputter coater (Quorum Technologies, Ltd, Kent, UK) and imaged 
using a Hitachi S3000N variable pressure SEM (Hitachi, Troy, MI). For histology, approx- 
imately 15 infected copepods were fixed in Alcohol-Formalin-Acetic Acid (AFA), embed- 
ded in mass in paraffin, sectioned at 6-8 um, and stained with hematoxylin and eosin. 
The morphology of adult worms was examined using a compound microscope by plac- 
ing worms on a microscope slide and compressing slightly with a coverslip. This allowed 
accurate counting of loculi and testes. 

Specimens obtained from experimental trials and from wild caught round stingrays were 
fixed in 95% ethanol for molecular analyses. Total genomic DNA was extracted from these 
samples using the Qiagen DNeasy Kit (Qiagen, Valencia, CA) according to the manufac- 
turer’s instructions. A 572 base pair region of the cytochrome c oxidase I (COI) gene was 
amplified via the polymerase chain reaction (PCR) using Cestoda-specific primers nLCO 
(S'-TTT ACT YTR GAY CAT AAG CGT-3’), and BEN-5 (5’-AAG CAG AAC CAA 
TTT ACG ATC-3’), according to Reyda and Marques (2011). Successful PCR reactions, 
determined via electrophoresis, were cleaned using MultiScreen HTS plates (Millipore 
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Corporation, Bedford, MA) and sequenced via Laragen, Inc. (Los Angeles, CA). Align- 
ments to close relatives were performed in Sequencher v4.10.1 and distance matrices 
(Jukes-Cantor model) were created using Geneious v8.0.2. Sequences of R. urobatidium 
obtained from experimental infections were compared to R. urobatidium, Ruptobothrium 
ditesticulum (Appy and Dailey, 1977) Coleman et al., 2019 and R. gravidum Friggens and 
Duszynski, 2005 from wild caught U. halleri, collected during the current study. COI 
sequences generated during this study are available via GenBank (Acc. #MN416032- 
MN416047). 

Tapeworm ontogenetic and morphological terminology generally follows that of Chervy 
(2002) for larval stages (but see Discussion), (Caira et al. 1999) for scolex attachment or- 
gans, Conn and Swiderski (2008) for embryonic envelopes, Chervy (2009) for microtriches; 
Jensen and Bullard (2010) for larval types, and Hilliard (1960) and Bylund (1975) for hex- 
acanth hooks. Miller and Lea (1972) was used to identify fishes. 


Results 


In the 2013 survey of fish parasites in Anaheim Bay, 75.6% (N = 41) of round stingray 
were found to be infected with | to 53 adult R. urobatidium. The only other rhinebothri- 
idean species present included R. gravidum (85.4%) and R. ditesticulum (26.8%). Presump- 
tive R. urobatidium plerocercoids were found in the gall bladder and bile ducts of several 
species of small benthic fish including 22.7% (N = 22) of longjaw mudsuckers (Gillichthys 
mirabilis Cooper, 1864) 23.1% (N = 26) of arrow and shadow gobies and 6.3% (N = 32) of 
staghorn sculpins (Leptocottus armatus Girard, 1854). These larvae were absent in topsmelt 
[Atherinops affinis (Ayres, 1860)] (N = 30), California killifish (Fundulus paryipinnis Gi- 
rard, 1854) (N = 37), shiner surfperch (Cymatogaster aggregata Gibbons, 1854) (N = 16), 
yellowfin goby [Acanthogobius flavimanus (TYemminck and Schlegel, 1845)] (N = 9), and 
barred sand bass (Paralabrax nebulifer (Girard, 1854)) (N = 60). Rhinebothrium urobatid- 
ium was particularly common in the tidal ponds where in 2013, 100% of arrow goby were 
infected with 2 to 54 worms. In the fall of 2018, 93% of 43 arrow and cheekspot gobies 
collected in Anaheim Bay tidal ponds were found to be infected with 1 to 57 plerocercoids 
per fish. 

Apolytic gravid proglottids of R. urobatidium (Fig. 1A) in the spiral intestine of round 
stingrays were initially distinguished/segregated from other gravid proglottids by a com- 
bination of moderate size relative to the gravid proglottids of other tapeworm species, 
uniform white appearance, a smooth surface and an indistinct midbody genital pore. Ident- 
fication as R. urobatidium was confirmed by examining egg morphology using a compound 
microscope. When placed in seawater, eggs extruded from an opening midway along the 
proglottid (Fig. 1A). Eggs of R. urobatidium taken directly from gravid proglottids were 
nonoperculate and enclosed by a transparent outer coat with a robust gradually taper- 
ing terminally rounded filament extending from one pole (Figs. 1B, 2A). Within this clear 
membrane the hexacanth embryo was surrounded by inner and outer envelopes. Expelled 
eggs held in seawater did not attach to one another or become tanned. Embryos were ac- 
tive within the envelopes. Hexacanth hooks were not uniform in size or shape (Table 1, 
Fig. 2B). The medial (third) pair were longer than the lateral (first) or medio-lateral (sec- 
ond) pairs and had a more recurved blade tip, and the mediolateral pair were shorter than 
the medial pair and had a wider handle and blade. The lateral pair was the most delicate 
with a gently curved blade tip and thin handle. 
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Fig. 1. Photographs of the development of R. urobatidium in T. californicus, held at 21°C. (A) Gravid 
proglottid expelling eggs. (B) Egg containing hexacanth embryo under coverslip pressure. (C) Hexacanth 
embryo in copepod intestinal lumen and in hemocoel at 8 hours post-exposure. (D) Hexacanth embryo 
within an envelope in copepod hemocoel. (E) Procercoids within a membrane in the copepod hemocoel 
at 10 d post-exposure (dpe). (F) Proceroid within a membrane at 10 dpe. (G) Procercoids isolated from a 
copepod at 14 dpe. (H) Mature procercoid within a membrate at 14 dpe. Abbreviations: as — apical sucker, 
c — cercomer, cc — calcareous corpuscles, h — hexacanth embryo, ie — inner envelope, m—membrane, oe — 
outer envelope, pr — procercoid. 
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Line drawings of R. urobatidium egg, and heat-killed procercoids from T: californicus and pro- 
cercoids and plerocercoids from the arrow goby at 21°C. (A) Egg. (B) Hexacanth hooks. (C) Procercoid 
with cercomer at 10 d post-exposure (dpe). (D) Procercoid from arrow goby at 8 dpe. (E) Stages in both- 
ridial development of the plerocercoid. (F) Plerocercoid with bothridia partially retracted. (G) Habitus of 
plerocercoid at 60 dpe. Abbreviations: as - apical sucker, b - bothridia, ba - bothridial aperture, c - cercomer, 


cc - calcareous corpuscles, g - gland cells and ducts, e - excretory system, je - inner envelope, | - lateral hook, 
m - medial hook, ml - mediolateral hook, oe - outer envelope. 
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Table 1. Morphometrics in micrometers (j1m)* of eggs and hexacanth larvae of R. urobatidium. 
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Number 10 
| 144+ 7.0 
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Filament length ne 
lament leng 79-103 
= hee 14+ 1.0 
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edlal Fair 8.6—9.5 


*Mean + standard deviation over the range. 


Experimental transmission of R. urobatidium to copepods was conducted in excess of 20 
times, and in all cases 100 percent of copepods became infected with up to 48 larvae per 
copepod. Mortality of copepods was common in the first days of infection. Two attempts 
to infect ostracods were not successful. At eight hours post-exposure, hexacanth embryos 
were present in copepod tissues lining the digestive tract, where they were surrounded by 
a clear space and thin membrane (Figs. 1C, D). Larvae penetrated the gut along its entire 
length but were found less frequently in the urosome (hindgut). As larvae grew they elon- 
gated, and by 5 dpe they had started to develop a cercomer and an apical organ/sucker. 
By 8 dpe calcareous corpuscles were present, the apical sucker appeared fully developed, 
and the cercomer was oval and connected to the body proper (Figs. 1F, G, H; 3A, B) by a 
narrow isthmus (Fig. 3C). Procercoids were clustered around the copepod digestive tract 
and surrounded by a thin membrane (Figs. 1E, F, H; 3B, C). Hexacanth hooks were not 
visible on the cercomer, but were occasionally observed in the posterior body of mature 
procercoids. While growth of some larvae continued through at least 20 dpe (Fig. 4A), 
elongation slowed following 8 dpe and by 11 dpe some larvae were near the maximum size 
of larvae examined at 20 dpe. Glandular cells were present in the anterior midbody of the 
procercoid, with ducts extending into the area of the apical sucker (Fig. 2C). The excretory 
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Fig. 3. Scanning electron micrographs of R. urobatidium in experimental (A-E, G-H, J—K) and natu- 
ral infections (F, I). A. Procercoid from T. californicus at 15 d post-exposure (dpe). B, C. Procercoid from 
T. californicus with enclosing membrane at 15 dpe. D. Plerocercoid from arrow goby at 30 dpe. E, G. Plero- 
cercoid from arrow goby at 51 dpe. F, H. Plerocercoid from natural infection of arrow goby. I. Adult scolex 
of specimens from round stingray. J. Bodies of plerocercoids extending from bile duct into the intestinal 
lumen of arrow goby. K. Procercoid microtriches at 15 dpe. L. Plerocercoid microtriches at 51 dpe. Abbre- 
viations: as - apical sucker, b - bothridium, ba - bothridial aperture, c - cercomer, ca - cercomer attachment, 
f - filithriches, 1 - intestine, m - membrane, s - spinitriches. 
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Fig. 4. Growth of R. wrobatidium in intermediate hosts held at 21°C: A. Procercoids in 7? californicus. 
B. Procercoids/plerocercoids in arrow goby. Solid line is total length of larvae while dashed line is change 
in length of bothridia. 


system consisted of a posterior common duct, which bifurcated in two lateral branches 
extending into the anterior end of the larvae. During development in the copepod, larvae 
usually did not move freely within the hemocoel, but adhered to the external surface of 
the digestive tract (Fig. 1E). No obvious bothridial development occurred within 20 dpe. 
However, one of two larvae present in a copepod at 30 dpe had developed rudimentary 
bothridia and was moving in the hemocoel. 

Gobies experimentally infected with R. urobatidium procercoids from copepods har- 
bored 11 to 44 larvae. At 5 dpe, procercoids with the apical sucker, but without a cercomer 
(Fig. 2D), were present between villi in the anterior intestine of experimentally infected go- 
bies (Fig. 5A). By 7 dpe some larvae, which had increased in length (Fig. 4B), were present 
in the bile duct and gall bladder (Fig. 5B); those in the gall bladder had attached to its 
inner wall by the apical sucker. Bothridia began forming in some larvae by 10 dpe as indis- 
tinct anterior swellings (Fig. 5C), and subsequently became septate oval bothridia, which 
developed anterior and posterior lobes as they increased in size (Figs. 2E; 3D, E, G; 5D). 
Septa at the anterior end formed first, with septa increasing in number and definition as the 
infection progressed (Fig. 2E). By 51 to 60 dpe, 36 to 38 loculi were present in the longer 
plerocercoids (Figs. 2E, G; 3G; 5E). During development, the anterior bothridial lobe was 
larger than the posterior lobe, and more loculi were present in the anterior lobe than in 
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Fig. 5. Photographs of R. urobatidium larvae in arrow goby held at 21°C: (A) Procercoids among the 
goby intestinal villi. (B) Procercoids attached to the inner wall of the gall bladder at 10 d post-exposure 
(dpe). (C) Procercoid beginning to differentiate into a plerocercoid at 15 dpe. (D) Plerocercoid with devel- 
oping bothridia at 30 dpe. (E) Developed bothrium attached to coverslip at 60 dpe. (F) Digestive tract of 
goby with plerocercoids in the gall bladder and bile duct. (G) Plerocercoid with bothridia and apical sucker 
retracted at 60 dpi. (H) Aperture of a retracted bothridium. I. Extended plerocercoid from a naturally in- 
fected arrow goby. Abbreviations: as - apical sucker, b - bothridia, ba - bothridial aperture, bd - bile duct, cc - 
calcareous corpuscles, g - gall bladder, i - intestine, im - inner envelope, oe - outer envelope, m - microtriches, 
pr - procercoid, pl - plerocercoid, s - stomach. 
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Table 2. Morphometrics in micrometers (jm)* of larval stages of R. urobatidium in copepods 
(7. californicus) and arrow gobies. ; 


Procercoid Plerocercoid Plerocercoid 
Parameter (20 dpe)! (51-60 dpe)! Natural Infection 
Number 14 6 7 
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Apical sucker width 36.44 wane AqD8 
Cercomer length at 20:0 s E 
21-46 
362+ 4.4 
Cercomer width ay = * 
Bothridium length Bese Nh oz 200.6 
170-387 213-382 
ig é 130 + 45.3 Ol se 754 
Bothridium width 98-900 ean Oni 
Number of loculi 32-38 34-40 


* Mean + standard deviation over the range. 
' dpe = days post-exposure. 


the posterior lobe (Figs. 2E; 3E). In heat-killed specimens, the shape and reflexion of the 
scolex/bothridia was similar to that in heat-killed plerocercoids and adult worms from 
natural infections (Figs. 3G, I). By 51 to 60 dpe, some plerocercoids approached the length 
of larvae present in naturally infected gobies (Table 2). Plerocercoids were densely packed 
with calcareous corpuscles and had excretory canals similar to those seen in the procer- 
coid. In heavy experimental and natural infections, the gall bladder and the bile duct were 
congested with plerocercoids (Fig. 5F) and in one case, the posterior portion of the larvae 
extended out of the bile duct into the intestine (Fig. 3J). In some cases, procercoids (lar- 
vae without bothridial development) were still present in the bile duct and gall bladder as 
late as 30 dpe. The bodies of plerocercoids were highly elastic (Fig. 51), and the bothridia 
and apical sucker could individually be retracted into the scolex (Figs. 2F; 3H; 5G, H). 
The pedicel retracted into the scolex followed by the posterior portion of bothridium and 
then the anterior portion. Plerocercoids viewed in the gall bladder of fish invariably were 
attached to the inner bladder wall with the apical sucker and with bothridia retracted. In 
both experimental and natural infections, developed plerocercoids were present in the bil- 
lary system but not in the intestine. No plerocercoids that could be attributed to a different 
morphotype were found in the bile/cystic duct or gall bladder of gobies although plero- 
cercoids of other tapeworms (believed to be Acanthobothrium spp. and Parachristianella 
Sp.) were present in the lumen of the posterior intestine and rectum. Mosquitofish did not 
become infected in two trial exposures. 

Presumptive filitriches were present on procercoids (Fig. 3K) but absent on the cercomer. 
In plerocercoids, filitriches were present adjacent to the apical sucker, and spinitriches were 
present on the corpus sometimes mixed with filitriches (Fig. 3L). Microtriches became 
Sparse toward the posterior of the plerocercoid. 
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Table 3. Genetic relationship of R. urobatidium (R.u.) from experimental infections with natural in- 
fections of R. urobatidium, at different localities, R. gravidum (R.g.) and R. ditesticulum (Ru.d.) based on 
mitochondrial cytochrome oxidase |. Numbers represent the difference in base pairs in a 465 base pair 
region. 


R. urobatidium 


ae ——ee——eE——EeEeEE—EE—E—E—————————e 
adult adult adult adult adult plero. plero. _ plero. 
Species Stage N* Locality’ AB LAH  LBH SB SDB AB 30dpe 60 dpe 
Ru.d. adult 1 AB — 123. 117-118 117-118 117-118 117 117 117 
R. g. adult 1 LAH - 117 117-119 117-118 118 117 NY 
R.u. adult 4 LBH 0-2 |-3 0-1 1-2 0-1 0-1 
adult 2 SB 3 1-3 1-2 1-2 1-2 
adult 4 SDB 0-2 ]—2 0-1 0-1 
plero. 1 AB — ] 
ER, 2 30 0 0 
plero. 1 60 = 


* N = number of specimens sequenced. 
' AB = Anaheim Bay, LAH = Los Angeles Harbor, LBH = Long Beach Harbor, SB = Seal Beach, 
SDB = San Diego Bay; plero. = plerocercoid. dpe = days post exposure. 


Cytochrome c oxidase I sequences of experimental plerocercoids were nearly identical to 
those from plerocercoids and adult R. urobatidium from natural infections, with only 1-3 
basepair differences out of a total 465-bp region, compared to 117-119 bp mismatch with 
other rhinebothriidean species (Table 3). 


Discussion 


The three-host life cycle of R. urobatidium includes a first intermediate host copepod 
where a caudate procercoid with an apical organ/sucker and cercomer develops, a second 
intermediate host small benthic fish in which a nonlacunate plerocerocoid bearing indi- 
vidually retractable bothridia and apical sucker develops in the biliary system, and the 
final host, the round stingray. Eggs of R. urobatidium include inner and outer envelopes, 
which are likely homologous with embryonic envelopes recognized by Conn and Swiderski 
(2008). The outermost clear membrane with a single rounded process may be the outer 
coat formed in some cestodes (Conn and Swiderski 2008). Other rhinebothriidean species 
including R. abaiensis Healy, 2006 (Healy 2006) and R. gravidum and R. ditesticulum (pers. 
obs.) have a transparent outer coat, although the eggs of these species all possess a process 
extending from each end of the egg. Other tapeworm species have similar outer coats (e.g. 
Proteocephalus spp. of the order Onchoproteodephalidea (Scholz 1999). The size and shape 
of these eggs is certainly a microecological adaptation that enhances transmission (Jarecka 
1961) either by providing buoyancy, being attractive/refractive, and/or by mimicking a 
particular food (e.g. phytoplankton) of the first intermediate host. 

The development of a procercoid with a cercomer and apical organ/sucker is a com- 
mon ontogenetic feature among the Cestoda (Freeman 1973, Jarecka 1975, Chervy 2002), 
but within the Rhinebothriidea this stage has only been described for Pseudanthobothrium 
Baer, 1956 (Jarecka and Burt 1984). Unlike in R. urobatidium, the cercomer of Pseudan- 
thobothrium spp. contains large refractile bodies. In addition, the body (protoscolex) of 
Pseudathobothrium is enclosed with a clear membrane within which the larva moves freely 
and separately from the cercomer (Jarecka and Burt 1984). In contrast, a clear membrane 
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observed in R. urobatidum surrounds the entire procercoid, including the cercomer, and the 
body of the procercoid is directly attached to the cercomer. In R. urobatidum the movement 
of the procercoid is constrained within the membrane, perhaps by virtue of its attachment 
to the cercomer. 

Among acetabulate fish tapeworms, development of a procercoid with an apical sucker 
and a cercomer, as in R. urobatidium, 1s also found in Acanthobothrium olseni (Mudry and 
Dailey 1971), A. hispidum (Riser 1956), and similar to initial development in species of 
Proteocephalus Wineland, 1858 (Scholz 1999). In most Proteocephalus spp. the formation 
of cercomer and apical sucker is quickly followed by the development of four suckers sim- 
ilar to what is found in the adult worms (Scholz 1999). In contrast, R. urobatidium did not 
develop bothridia in the first 20 d in the copepod, and larvae without bothridia fed to gob- 
ies at 15 dpe, developed bothridia similar to the adult worm. Thus in R. urobatidium reared 
in T. californicus, the development of bothridia is not required in order to infect the sec- 
ond intermediate host as it apparently is in Proteocephalus. The pattern in R. urobatidum 
is more similar to that in the bothriate cestodes such as some bothriocephalideans or di- 
phyllobothriideans, in which the procercoid does not develop an adult-like scolex until it 
enters the second intermediate fish host (Freeman 1973). 

The presence of microtriches on procercoids has been observed in other cestode orders 
(Davydov et al. 1995). While not studied in detail herein, a transition of microthrix type 
from procercoid to plerocercoid of R. urobatidium has been reported for other tapeworms 
(Hopkins and Charles 1969). Such a transition is not surprising in light of general ontoge- 
netic tegumental changes that occur in other tapeworm species (Davydov et al. 1995). 

The membrane surrounding the procercoid of R. urobatidium appears to be of parasite 
rather than host origin and similar membranes have been observed in other cestode species. 
Invertebrate immune responses typically involve a cellular response and encapsulation (van 
der Veen and Kurtz 2002; Rowley and Powell 2007). Such responses were not evident in 
T. californicus infected with R. urobatidium. In contrast, ultrastructural and in vitro stud- 
ies have variably described a surface coat in the spathebothriidean Cyathocephalus trun- 
cates (Pallas 1781) (Okada 1990), an electron-translucent fibrillar secretion around the 
procercoid of the bothriocephalid Triaenophorus nodulosus (Pallas, 1781) (Davydov et al. 
1995), a surface filamentous coat in the caryophyllidea Archigetes sieboldi Leuckart, 1878 
(Poddubnaya et al. 2003), and an outer surface or outer layer in the diphyllobothriidean 
Schistocephalus solidus (Muller 1776) (Jakobsen et al. 2012; Marwaha et al. 2013). The 
membrane/layer sometimes described as carbohydrate (Jacobson and Doyle 1996) or shed 
outer surface (Jakobsen et al. 2012) has been purported to help the parasite evade the im- 
mune responses of the invertebrate intermediate host (Hammerschmidt and Kurtz 2005) 
or to provide protection from enzymatic action when the procercoid is passing through 
the vertebrate host digestive tract (Kurtz and Franz 2003; Marwaha et al. 2013). An ultra- 
structural/cytochemical study would be needed to better compare the membrane found in 
R. urobatidium with those found surrounding the procercoids of these other tapeworms. 
Absence of bothridial development in R. urobatidium in the copepod may also be related 
to the presence of this membrane. This is supported by our observation of a single larva of 
R. urobatidium with bothridia moving freely in the copepod hemocoel at 30 dpe. In addi- 
tion, at least rudimentary bothridia are present in larvae of R. gravidum and A. parvinunci- 
natum, which move freely in the hemocoel of 7: californicus (pers. obs.). This may suggest 
that conditions in the copepod hemocoel trigger or allow development of bothridia. 

The presence or absence of bothridia in the larvae in the first intermediate host may also 
be related to the trophic demands of transmission. In Proteocephalus spp., development 
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of acetabula/scolex in the first intermediate host precedes subsequent development to an 
adult in the digestive tract of the final host in a usual two-host life cycle (Scholz 1999). 
In acetabulate fish tapeworms in a three-host cycle, such as R. urobatidium, the absence of 
bothridial development may be related to subsequent location of the larva in the second in- 
termediate host; e.g., bothridia may be an impediment to larval parasites that migrate into 
ducts (e.g. R. urobatidium) or the body cavity [e.g., Type 5 larva of Chambers et al. (2000) 
and R. ditesticulum (pers. obs.)] of the second intermediate host. In any case, while the 
ontogenetic stages (procercoid and plerocercoid/merocercoid) of acetabulate tapeworms 
appear similar across taxa, the timing of larval morphogenesis from one to the other in 
the first intermediate host appears to be variable and likely a species-specific adaptation to 
trophic requirements of transmission. 

The plerocercoid of R. urobatidium conforms to Jensen and Bullard’s (2010) Type VII, 
which has an apical sucker and four stalked bothridia that are individually and totally 
retractable into the scolex. This general type of larva has been reported from the diges- 
tive system of a number of fish species from the Great Barrier Reef (Chambers 2000) and 
in the Gulf of Mexico (Jensen and Bullard 2010). The larval worms described by Brooks 
and Brothers (1997) from arrow goby and the shadow goby, Quietula y-cauda (Jenkins 
and Evermann, 1889) (Gobiidae) in San Diego Bay, California are most likely R. uroba- 
tidium. Small benthic fish in Anaheim Bay are commonly infected with R. urobatidium, 
especially in the man-made tidal ponds. The heavy use of these tidal ponds by gestating 
round stingray in the late summer and fall (Jirik and Lowe 2012) evidently provides a focus 
of infection of small benthic fish with R. urobatidium. The occurrence of rhinebothriidean 
larvae in fish is not always consistent with the feeding habits of final hosts, which consume 
few fish (Jensen and Bullard 2010). This is also true for R. urobatidium, in which the round 
stingray’s diet consists of only 7.5% fish (Babel 1967). 

While the presence of larval shark tapeworms in the gall bladder and bile/cystic duct 
of a fish host does not appear to be a common occurrence (e.g. see Jensen and Bullard 
2010) larvae that might be attributed to Rhinebothrium or related taxa have been reported 
from the biliary system of fish (Linton 1897, 1905; Chandler 1935; Overstreet 1968; Jensen 
2009). In tapeworms of groups other than elasmobranchs, acystic larvae of Valipora spp. 
(Cyclophyllidea: Gryporhynchidae), parasites of birds, are also found in the gall bladder of 
bony fishes (Jarecka 1970; Scholz and Salgado-Maldonado 2001). Larval helminths exploit 
intermediate host tissues to avoid the gut, maximize their growth, and minimize mortality 
(Chubb et al. 2009; Parker et al. 2009a, b). The specialized location of R. urobatidium in 
the biliary system, which might seem to be an inhospitable environment, would provide the 
benefits of being in host tissues or body cavities while avoiding the need to penetrate host 
tissues. In particular, R. urobatidium plerocercoids in the biliary system avoid energy expen- 
diture associated with withstanding gut peristalsis, possibly avoid host immune response 
due to minimal host/parasite contact, and evidently find adequate nutrition to significantly 
increase their size. The protruding of the posterior ends of larvae into the intestinal lumen 
could be a means to obtain nutrition (Fig. 3J) although this was only observed in heavy 
experimental infections of the arrow goby. 

The larvae of R. urobatidium, like those of at least some other Rhinebothrium species 
(see Jensen and Bullard 2010), possess neither evaginated bothridia (plerocercoid) nor an 
invaginated scolex (merocercoid). Instead, bothridia of R. urobatidium are individually 
retractable; they are withdrawn, without invagination of the bothridia, into the scolex. 
In discussing larval cestode terminology Chervy (2002) recognized the difficulty in the 
use of scolex retraction and invagination as a feature in defining larval types but felt it 
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may be a useful feature. Other authors have also noted difficulty in assigning a designa- 
tion to the bothridia of plerocercoids (Rocka 2003; Guagliardo et al. 2009). While larval 
R. urobatidium from gobies do not appear to fit neatly into the definition of a plerocercoid 
or merocercoid it seems premature to create new or resurrect terminology in light of the 
absence of transmission and ultrastructural studies on the thousands of shark tapeworm 
metacestodes. Therefore, we have utilized the term plerocercoid, which is consistent with 
the terminology used by Jensen and Bullard (2010) for similar larvae. 

The number of bothridial loculi present in the plerocercoid of R. urobatidium approaches 
or reaches the number of loculi in adult worms (Table 2; Appy and Dailey 1977), suggest- 
ing that in regions where the local elasmobranch tapeworms are known, plerocercoids of 
Rhinebothrium may be attributed to an adult form and thus a species based on the both- 
ridial morphology of the plerocercoid. In addition to morphological similarity between 
experimental plerocercoids and adult worms, and plerocercoids from natural infections, 
mitochondrial DNA sequence data confirmed the finding that among rhinebothriidean 
tapeworms commonly found as adults in round stingray (1.e., R. gravidum and R. dites- 
ticulum) in Anaheim Bay and Southern California, larvae from experimental infections 
are those of R. urobatidium (Table 3). Similarly Laskowski and Rocka (2014) found fewer 
than five base pair differences in COI] sequence between plerocercoids from marbled rock- 
cod, Notothenia rossii Richardson, 1844, and the onchoproteodephalidean Oncobothrium 
antarticum Wojciechowska, 1990 from the skate Bathyraja eatonii (Gunther 1876), a find- 
ing also confirmed with IsrDNA, which not surprisingly showed no molecular differences. 
In contrast, Jensen and Bullard (2010) found that congeners of some shark tapeworms 
may have identical COI sequences and recommended caution in identification of larvae 
or adults to species using COI. However, even absent genetic data, R. urobatidium can be 
distinguished from other rhinebothriidean tapeworms occurring in round stingray. Rup- 
tobothrium ditesticulum, which 1s relatively scarce in round stingray in Anaheim Bay, has 
bothridia comprising two separate lobes, very large gravid proglottids and eggs with fil- 
aments at both poles. The bothridia of R. gravidum are not indented/lobed, their gravid 
proglottids have crenulated surface, prominent genital atrium and eggs have filaments at 
both poles. Rhinebothrium urobatidium exhibits substantial growth in both intermediate 
hosts. Procercoid growth was asymptotic and reached growth arrest at larval maturity 
(GALM) (Ball et al. 2008; Parker et al. 2009a) between 8 and 11 dpe at 21°C (Fig. 4). 
In contrast, while the morphological development and bothridial growth of the plerocer- 
coids proceeded sequentially, growth in length of the larvae was quite variable. This is likely 
a fixation artifact related to highly contractile plerocercoid, and a less-than-optimum ex- 
perimental design including the examination of only a single goby at each sampling point, 
and possible over-infection of small (~2 cm) goby hosts. In particular, extremely high in- 
tensity of infection has been shown to cause delayed development in the intermediate host 
(Benesh 2010). The goby examined at 30 dpe in the present study (Fig. 4) was infected 
with 44 larvae, which completely filled the gall bladder and greatly expanded the width 
of the bile duct. It should also be noted that experimental plerocercoids at 51 to 60 dpe 
did not reach the maximum size of larvae from natural infections (Table 2). Additional ex- 
perimental studies will be necessary to better determine when growth arrest occurs in the 
plerocercoid and at what point in development plerocercoids become infective to the final 
host. 

While recent years have seen a rise in morphological/taxonomic and phylogenetic in- 
formation on elasmobranch tapeworms (e.g. Caira and Jensen 2014; Caira et al. 2014; 
Caira and Jensen 2017), there are few experimental/ultrastructural studies on any of the 
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nine orders of elasmobranch tapeworms. This is in contrast to the significant advances in 
experimental and ultrastructural studies of some mammalian and avian tapeworms (e.g. 
Michaud et al. 2006; Hammerschmidt and Kurtz 2009) including the sequential in vitro 
development of the procercoid and pleroceroid of Schistocephalus solidus (Jakobson et al. 
2012). The present study experimentally confirms the identity of elasmobranch tapeworm 
plerocercoids found in small benthic fish as R. urobatidium, and for the first time describes 
morphogenesis of larval stages from first to second intermediate hosts. However, it still 
falls short of documenting the exact transmission that occurs in nature. It seems clear that 
R. urobatidium has a three-host life cycle, but the species of the natural first invertebrate 
(copepod) intermediate host(s) of R. urobatidium is still unknown, and worms have yet to 
be experimentally transmitted to the final host to determine time to maturity and longevity. 
Finding copepods infected with R. urobatidium in the marine environment is a difficult un- 
dertaking (although not impossible — see Marcogliese 1995), and use of wild-caught round 
stingrays for experimental studies is not feasible since regardless of size or locality of cap- 
ture, they are invariably infected with a menagerie of tapeworm species. We are currently 
studying the use of oral anthelminthics to remove digestive tract parasites from wild-caught 
stingrays and husbandry of stingrays born live in captivity to allow more complete trans- 
mission studies. Round stingray is a particularly good host for experimental studies since 
it is well-studied locally (Babel 1967; Hoisington and Lowe 2005; Vaudo and Lowe 2006; 
Mull et al. 2008; Plank et al. 2010; Jirik and Lowe 2012; Lyons et al. 2017), is abundant 
and easily caught in shallow waters where large congregations of rays result in foci of para- 
site transmission, is small and relatively easy to maintain/feed in captivity, and above all, it 
has a diverse tapeworm fauna including representatives from a number of tapeworm orders 
(Rhinebothriidea, Trypanorhyncha, Diphyllidea, and Onchoproteocephalidea). 
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Abstract.—Seagrasses are an important refuge for fishes and provide ecosystem ser- 
vices worldwide. Along the Pacific Coast, however, quantitative assessments of the 
ecological role of seagrass habitats with limited anthropogenic impacts are lacking de- 
spite their value to resource management. To address these issues, the ecological func- 
tion of the eelgrass Zostera marina Linnaeus in a protected area off Santa Catalina 
Island, California, USA was quantified using estimates of secondary production in 
the kelp bass Paralabrax clathratus (Girard, 1854). Monthly assessments of Zostera 
structural complexity, as well as the size and abundance of juvenile kelp bass were 
used to establish a baseline of fish biomass and recruitment associated with Zostera 
habitat. The greatest number of kelp bass was recorded in the summer and fall months 
and the fewest during the winter and spring. Secondary production in kelp bass (0.1 to 
0.59 gm * mo!) followed monthly changes in Zostera habitat structure throughout 
the 24-month study period. Seasonal changes in the structural complexity of coastal 
Zostera beds may influence the flow of energy to adjacent reef environments through 
the transfer of fish biomass. These findings help define the potential role of Zostera 
habitat within the larger context of nearshore coastal ecosystems in the Southern Cal- 
ifornia Bight. 


Coastal seagrass beds are an important nursery habitat for fishes and provide a variety 
of ecosystem services worldwide (Worthington et al. 1992; Heck et al. 2003; Burkholder 
et al. 2007; Warren et al. 2010; Parsons et al. 2015). Given the significance of this resource, 
the protection of seagrass habitat has received considerable attention (Coles et al. 2014; 
Bas Ventin et al. 2015; Schultz et al. 2015). The National Oceanic and Atmospheric Ad- 
ministration’s National Marine Fisheries Service (NOAA Fisheries) recognizes seagrasses 
as habitat areas of particular concern, which are defined as subsets of essential fish habitat 
under the Pacific Coast Groundfish Fishery Management Plan. Such areas are considered 
high priorities for conservation and management, and warrant special attention during the 
regulatory process set forth by the Magnuson Stevens Fishery Conservation and Manage- 
ment Act in 1976. Nevertheless, the management of seagrasses has proven challenging in 
some areas, particularly those with remote coastlines that are difficult to access, and in 
other cases due to a paucity of information on seagrass resources altogether (Long and 
Thom 2001). For example, because of their isolation from the Southern California main- 
land, the distribution and habitat function of seagrass beds off the California Channel 
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Islands are poorly known; this is in contrast to the many long-term studies of shallow 
rocky reefs and kelp forests from the same region (Engle and Miller 2005). 

One of the most common species of seagrass found along the Southern California coast 
is the eelgrass Zostera marina Linnaeus 1753, which forms extensive beds in shallow bays 
and estuaries with low to moderate current exposure (Dawson and Foster 1982; Engle 
and Miller 2005). Many of the areas inhabited by Zostera are subject to anthropogenic 
stressors such as dredging, eutrophication and development that decrease the ecological 
function and stability of these communities (Obaza et al. 2015; Shelton et al. 2017). While 
studies of Zostera beds found in these locations offer information on habitat function in 
a stressed environment, data from protected areas provide an upper bound of the behav- 
ioral response by fishes that are dependent on seagrass as a refuge from predators (Jackson 
et al. 2001; Gillanders et al. 2003). Such is the case for populations of Zostera around the 
California Channel Islands, which are separated from many of the anthropogenic stressors 
found along the mainland coast and thus are an ideal location to collect ecological and en- 
vironmental data (Findlay and Allen 2002; Davis 2005; Mason and Lowe 2010; Saarman 
and Carr 2013). Of particular interest is information on habitat use and movement patterns 
of ecologically and economically important populations of reef fishes living in these areas 
(NOAA 2014). Accordingly, estimates of secondary production are used to quantify the 
biological benefit of different habitats by accounting for both changes in fish biomass and 
population parameters that include fish mortality, recruitment, immigration and emigra- 
tion (DeMartini et al. 1994; Johnson et al. 1994; Love et al. 1996; Faunce and Serafy 2008; 
Jeong et al. 2009; Kamimura et al. 2011; Yeager et al. 2012; Williams et al. 2013; Claisse 
et al. 2014). 

Here, we compare secondary production in juvenile kelp bass Paralabrax clathratus 
(Girard 1854), a temperate serranid fish, in relation to the three-dimensional (= struc- 
tural) complexity of Zostera habitat in a protected area off Santa Catalina Island, Cal- 
ifornia, USA. Zostera forms highly productive beds (Risgaard-Petersen et al. 1998) that 
function as an important nursery habitat for kelp bass (Mendoza-Carranza and Rosales- 
Casian 2002; Altstatt et al. 2014), as well as other common and economically important 
fishes such as rockfish (Sebastes sp.) and surfperch (Embiotocidae) (Hoffman 1986; Allen 
et al. 2002). Our findings provide an important baseline of fish biomass and recruitment 
associated with Zostera habitat in the absence of significant anthropogenic impacts. By 
integrating multiple metrics for evaluating Zostera habitat function, such as species use 
and secondary production, we aimed to improve conservation and management efforts by 
evaluating the role of Zostera habitat in a Southern California marine protected area. 


Material and Methods 
Study Site and Natural History 


Field surveys were conducted monthly from September 2013 to August 2015 in Big Fish- 
erman’s Cove (33°26'38.85’N, 118°29'6.86”W) in the Blue Cavern State Marine Conserva- 
tion Area (SMCA) on Santa Catalina Island, ~35 km off the coast of Southern California, 
USA (Fig. 1). Neither anchoring nor fishing is allowed in Blue Cavern SMCA, and all com- 
mercial and recreational take of marine life is prohibited as outlined by California’s Marine 
Life Protection Act of 1999. The subtidal habitat of Big Fisherman’s Cove extends seaward 
from a northwest facing, boulder and rock cobble beach as a sloping, sandy bottom that 
steadily descends to depths >35 m that lie towards the center and mouth of the cove. The 
northeast and southwest borders are comprised of semi-vertical rocky reef (5-20 m depth) 
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Fig. 1. Map showing the location of (a) Santa Catalina Island relative to the Southern California main- 
land. Surveys of (b) eelgrass Z. marina and kelp bass P. clathratus were conducted in (c) Big Fisherman’s 
Cove, which is part of the Blue Cavern State Marine Conservation Area. Hatched area indicates survey site 
and extent of Zostera areal coverage (~0.24 hectares). Map adapted from Froeschke et al. (2006). 


covered with macroalgae and kelp that provides a forest habitat for fishes and invertebrates 
(Abbott and Hollenberg 1992; Parnell et al. 2010). The sandy, inner reaches of the cove 
(4-12 m depth) support a Zostera bed that covers an area ~0.24 ha in size (Obaza and 
Ginsburg unpublished data). Areal coverage of Zostera in Big Fisherman’s Cove was esti- 
mated using a Trimble® R1 Global Navigational Satellite System (GNSS) receiver held at 
the water surface in which a scuba diver was tracked while swimming along the perimeter 
of the bed at depth. Patches of Zostera were first documented at this site in 1996 and are 
likely the result of natural colonization from nearby populations (Engle and Miller 2005). 

Fish surveys targeted kelp bass, which are one of the most important recreational species 
of reef fish off Southern California (Young 1963; Erisman et al. 2011). Typically found 
in shallow water (3-25 m depth), juveniles feed on benthic invertebrates and switch to a 
mostly piscivorous diet as adults (Hobson and Chess 1976; Eschmeyer and Herald 1999). 
Mature individuals form breeding aggregations in the late spring to early fall in which 
larvae enter the plankton and settle after 28-30 days onto Zostera beds (Valle et al. 1999; 
Allen et al. 2002) and shallow, rocky reefs (Carr 1994; Love et al. 1996; Cordes and Allen 
1997; Erisman and Allen 2006). Adult kelp bass exhibit high levels of site fidelity (Mason 
and Lowe 2010). However, catch and release studies by Carr (1994) and Hartney (1996) 
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indicate that juveniles are unlikely to return to the reef in which they settled (home) after 
moving to a new location. 


Zostera and Kelp Bass Surveys 


The structural complexity of Zostera habitat, as well as the size and abundance of juve- 
nile kelp bass was measured by scuba divers. Surveys were conducted along four separate, 
30-m long benthic transects spaced ~4 m apart. The beginning of each transect was marked 
with a semi-permanent sand anchor and surveyed along a northwest compass heading of 
310 degrees. Structural complexity was defined using three key measurements of Zostera 
habitat: shoot frequency, shoot density and canopy height. Frequency was estimated using 
a point-intercept approach in which Zostera leaf shoots were recorded as either present or 
absent along each meter interval (from | to 30 m) directly beneath the transect line. Density 
was measured as the total number of leaf shoots counted within four quadrats (0.25 m7) 
placed at 10-m intervals (0-30 m) along each transect. Lastly, Zostera canopy height was 
recorded (to the nearest cm) as 80% of the mean length of >10 haphazardly selected leaf 
shoots (Short and Durate 2001) within each of the quadrats described above. The product 
of each of these metrics (i.e., frequency, density and canopy height) was used to create an 
eelgrass structural index in order to view the collective changes in Zostera habitat over time. 

Juvenile kelp bass were surveyed visually from 1 m above the benthos and within a 2 m 
wide swath along each 30 m transect line. Given the importance of Zostera habitat to new 
recruits (Jackson et al. 2001), visual surveys focused on juvenile kelp bass (2.1—10 cm TL; 
Love et al. 1996) and not adults that may have migrated from other areas. The size and 
abundance of juveniles was recorded by divers as they swam the length of each transect 
line. The TL of an individual was recorded to the nearest cm following the methods of 
Bell et al. (1985). Kelp bass abundance was measured as the total number of juvenile fish 
recorded on a given transect. Visual estimates of fish length and density are known for 
a range of different reef fishes and habitats (Sale 1980; Brock 1982; Coyer and Witman 
1990; Ebeling and Hixon 1991) and are frequently used to monitor subtidal marine life off 
Southern California (Pondella et al. 2006; Gillett et al. 2012; Coates et al. 2018). 


Kelp Bass Secondary Production 


Secondary production was based on visual surveys of length and abundance of juvenile 
fish using a modified version of the model developed by Claisse et al. (2014). Specifically, 
secondary production rates were estimated for individual kelp bass observed on a given 
survey by determining changes in biomass after one month of growth using a weight— 
length relationship specific to this species. A linear growth function was used, rather than 
an exponential growth curve, as the former provides a better fit for modeling the growth of 
young-of-the-year (YOY) fishes (Faunce and Serafy 2008; Yeager et al. 2012; J.T. Claisse 
pers. comm. ). 

Fish biomass (i.e., mean wet body weight, W) was calculated from the observed length 
of kelp bass using the weight-length relationship: W = a TL”, where TL is the total length 
of individual kelp bass from visual surveys and the constants a and 5 (0.00813 and 3.03, 
respectively) are specific growth rates reported for this species (see Froese et al. 2014). 
The age (months-old) of individual kelp bass was determined using empirical age-at-length 
values for juveniles (2.1 to 10 cm TL) reported by Love et al. (1996). From these data, a 
linear growth rate of 1.17 cm TL mo7! was calculated for post-settlement (>2.1 cm TL) 
kelp bass <1 year old. Individual kelp bass, whose TL was recorded (to the nearest cm) 
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on dive surveys in Big Fisherman’s Cove, were then sorted into one of the following age- 
at-length bins (AL;—ALg): | month old, 2.1—3.27 cm TL; 2 months old, 3.28—4.45 cm TL; 
3 months old, 4.46—5.63 cm TL; 4 months old, 5.64—6.81 cm TL; 5 months old, 6.82—7.99 
cm TL; 6 months old, 8.0-9.17 cm TL). The growth (G) of individual kelp bass from each 
length bin was calculated using visual survey data as the difference between changes in 
biomass, estimated by TL, after one month of growth (ATL = 1.17 cm mo7!) using the 
equation: G = a(TL + ATL)? — W. 

The instantaneous rate of fish mortality (M) was calculated monthly with the TL of 
individual kelp bass from each age-at-length bin using the empirical formula reported by 
Gislason et al. (2010) as follows: In (M) = 0.55 — 1.61 [In (TL)] + 1.44 [In (Lo)] + In 
(k), where Lo and k (69.8 cm TL and 0.06, respectively) are key life-history parameters 
specific to kelp bass reported by Love et al. (1996). We selected this approach because it 
is the best-supported estimator for determining the number of individuals that survive to 
successive life stages (Lv and Pitchford 2007; Claisse et al. 2014). Fish mortality rates were 
estimated on a monthly basis (as opposed to an annual schedule). Survivorship (S) of kelp 
bass was then calculated using the monthly length- and species-specific exponential rate of 
fish mortality for each month sampled using the equation: S = e ™. 

Accumulated biomass (i.e. somatic production, P,) of kelp bass was calculated on a 
monthly basis as the product of fish growth (G), survivorship (S) within a given age-at- 
length bin and fish density (N; number of juvenile kelp bass recorded within a 2 m wide 
swath along each of the four, 30 m transects [ = 240 m? total survey area]) using the equa- 
tion: P, = (G) (S) (N). Total secondary production (Pr) was calculated as the monthly sum 
({) of recruitment (Pr, accounts for YOY fish growth) and somatic production (Ps) for 
each of the six post-settlement age-at-length bins (AL;-AL¢) recorded over the course of 
the study as follows: Pp = Xm (Pr + Pari + Par2 + Par3 + Para + Pars + Pato). 


Data Analysis 


The relationship between the structural complexity of Zostera habitat and the abun- 
dance of juvenile kelp bass was analyzed using a multiple regression model. Visual esti- 
mates of fish density were tested for heteroscedasticity and treated as individual, monthly 
events. Interactions among the different predictor variables of Zostera habitat (leaf shoot 
frequency, density and canopy height), measured across different months, were presumed 
to have at least one common variable between them. Temporal autocorrelation among 
covariate months and years was accounted for by applying a generalized least squares 
(GLS) and maximum likelihood (ML) estimation procedure to all statistical models tested. 
Akaike’s Information Criterion (AIC) was used to determine the model of best fit. Com- 
parisons between the different predictor variables were tested for multicollinearity using 
a Farrar chi-square test. Although the relationship between Zostera shoot frequency and 
canopy height was collinear (x7 = 2.39, p < 0.05), the relationship between each of these 
factors in addition to shoot density was not strong enough to warrant a statistical correc- 
tion. Furthermore, variance inflation factors for each of the Zostera predictor variables 
were not statistically significant (critical value < 2). Therefore, multicollinearity was un- 
likely to pose an issue (Graham 2003) for estimating the relationship between the structural 
complexity of Zostera and the presence of kelp bass. 

Comparisons between the abundance of kelp bass and each of the predictor vari- 
ables of Zostera habitat structure were analyzed separately using a GLS approach, and 
subsequently evaluated by AIC using the information-theoretic approach for analyzing 
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Fig. 2. Monthly changes in (a) frequency, (b) density and (c) canopy height of eelgrass Zostera marina 
recorded from September 2013 to August 2015 in Big Fisherman’s Cove. Values are means + | SE. Where 
not shown, errors fell into the graphical representation of the data point. 


ecological data developed by Burnham and Anderson (2002). Given the relatively low 
sample sizes, a second-order (i.e., corrected) information criterion (AICc) approach was 
used in place of the general AIC estimator. These values were, in turn, compared to kelp 
bass abundance and secondary production using a GLS approach. Error values reported 
are += | SE of the mean. All data were analyzed using R Statistical Software (R Core Team 
2015) with the dplyr (Wickham et al. 2015) and n/me (Pinheiro et al. 2015) packages. 


Results 


The frequency of Zostera leaf shoots (Fig. 2a) increased nearly 3-fold from March 
to July 2014 (average = 7.4 + 0.6 to 21.0 + 0.78 shoots m7!, respectively; ANOVA, 
F, 9 = 8.13, p < 0.05) and remained relatively unchanged until December 2014 (average = 
ANOVA, F}, 1; = 26.9, p = 0.0004). Shoot frequency then dropped 32% from an average of 
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Table 1. Factor type (Ws), AIlCc scores, AIlCc differences (A) and Akaike weights (w) for predicting the 
use of Zostera habitat structure by juvenile kelp bass Paralabrax clathratus. The model with the smallest 
AAICe value was selected as the best fit: ¢ are the maximized log likelihoods, Cumul (w) are the cumulative 
Akaike weights and Coeff (W) are the factor coefficients. 


W AlICc AAICc w Cumul () Coeff (Ws) 
Canopy height 257.48 0 0.37 -122.08 0.37 1.56 
Shoot density 0.84 
Canopy height YXSoO 1.06 0.21 -124.22 0.59 1.81 
Canopy height 259.45 OY 0.14 -123.06 0.73 1.43 
Shoot frequency 0.97 
Canopy height LIDS ool) Ome -121.46 0.84 33) 
Shoot density 0.72 
Shoot frequency 0.69 
Shoot density 260.96 3.47 0.07 -123.81 0.90 0.87 
Shoot frequency 1eS2 
Shoot density 261.31 309) 0.05 -123.60 0.96 0.32 
Shoot frequency 261.90 4.42 0.04 -125.90 1.00 1.81 


19.9 + 0.4 shoots m7! between July and December 2014 to 13.6 + 0.33 shoots m7! between 
January and August 2015 (ANOVA, F. \3 = 26.1, p = 0.0003). Alternatively, the density 
of Zostera leaf shoots (Fig. 2b) increased 39% from an average of 22.2 + 2.54 shoots m-? 
between September 2013 and June 2014 to 30.9 + 1.93 shoots m~’ between July 2014 to 
May 2015 (ANOVA, F290 = 7.7, p < 0.05). From here, shoot density more than tripled in 
December 2015, reaching a peak value of 110.3 + 24.1 shoots m~* (ANOVA, F;. 13 = 21.5, 
p = 0.0006). Finally, despite month-to-month fluctuations in the canopy height of Zostera 
(Fig. 2c), the overall length of leaf shoots did not significantly change throughout the study 
period (average = 20.0 + 1.13 cm; ANOVA, F 93 = 0.33, p = 0.57). 

A total of 1,419 kelp bass were recorded within the Zostera study site in Big Fisher- 
man’s Cove (Fig. 3). Nearly 87% of fish were identified as juveniles (2.1 to 10 cm TL; Love 
et al. 1996) with 84% of individual kelp bass estimated as younger than 6 months of age 
(<9.18 cm TL). Kelp bass were not observed in the Zostera study area on the following 
dates: February 2014, April 2014, January 2015 and May 2015. The number of kelp bass 
recorded for each age-at-length bin ranged from 124 to 311 individuals (average = 202.7 + 
26.8 kelp bass) over the course of the study (Fig. 3a-g). Kelp bass abundance across all age- 
length bins varied seasonally with the greatest number of fish recorded in the summer and 
fall months (June to November) and the fewest during the winter and spring (December to 
May) averaging 98.3 + 14.9 and 19.9 + 9.3 kelp bass mo™!, respectively (ANOVA, F\. 23 = 
19.9, p = 0.0002). 

Seasonal changes in Zostera habitat structure (Fig. 4a) were observed throughout the 
study period reaching its peak in the summer and fall (June to November, average = 
15,519 + 3,099) and dropping to its lowest point in the winter and spring (December 
to May, 6,381 + 1,838; ANOVA, F; 23 = 6.34, p < 0.05). Likewise, seasonal changes in 
Zostera habitat structure were significantly correlated with juvenile kelp bass abundance 
(BGtm = 0.015, p < 0.001; Fig. 4a). Results of the AICc goodness-of-fit model selection 
are shown in Table 1. The model that best approximated the use of Zostera habitat struc- 
ture by kelp bass, as indicated by the smallest AAICc score, included eelgrass shoot density 
and canopy height as the predictor variables (AICc = 257.5). Two other models accepted 
by the AIC, statistic (AAICc < 2) also included Zostera canopy height as a predictor 
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Fig. 3. Monthly abundance of kelp bass P. clathratus observed amongst eelgrass Z. marina habitat in 
Big Fisherman’s Cove. A total of 1,419 kelp bass were recorded from September 2013 to August 2015. Kelp 
bass observed during visual surveys were sorted into different age-at-length bins as follows: (a) 1 month 
old, 2.1-3.27 cm TL, (b) 2 month old, 3.28—4.45 cm TL, (c) 3 month old, 4.46—5.63 cm TL, (d) 4 month 
old, 5.64-6.81 cm TL, (e) 5 month old, 6.82-7.99 cm TL, (f) 6 month old, 8.0-9.17 cm TL and (g) > 9.18 
cm (subadult stage and older). On February and April 2014 and January and May 2015, no kelp bass were 
observed in the Zostera study area. 
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Fig. 4. Relationship between the (a) structural complexity of eelgrass Z. marina habitat and kelp bass 
P. clathratus abundance, and (b) rates of secondary production in juvenile kelp bass surveyed monthly from 
September 2013 to August 2015 in Big Fisherman’s Cove. Eelgrass structure index defined as the product of 
the mean leaf shoot frequency, density and canopy height of Zostera. 


variable: both on its own (AICc = 258.6) and when combined with leaf shoot frequency 
(AICc = 259.5). AAICc scores for the remaining candidate models exceeded the threshold 
value for the best-fit model (AAICc < 2) and, therefore, were not considered an accurate 
predictor of the relationship between kelp bass and Zostera habitat structure. 

Rates of secondary production calculated for juvenile kelp bass (range = 0.1 to 
0.59 g m-? mo7!) in Big Fisherman’s Cove were significantly correlated with monthly 
changes in Zostera habitat structure throughout the study (BGrm = 0.21, p = 0.0005; 
Fig. 4, Table 2). Secondary production in kelp bass was 6-times higher in the summer 
and fall months (June to November, average = 0.3 + 0.05 g m-° mo~!) than in the winter 
and spring (December to May, average = 0.05 + 0.02 g m-* mo7!). However, interannual 
shifts in fish secondary production from 2013 to 2015 were not observed (range = 1.62 + 
0.92 to 2.39 + 0.70 gm yr7!; Bom = 0.56, p = 0.07; Table 2). 


Discussion 


Data presented here provide an important baseline of fish biomass and recruitment as- 
sociated with Zostera habitat in the absence of significant anthropogenic impacts, and will 
help to define the role of seagrass within the larger context of nearshore coastal ecosys- 
tems. Because juvenile fishes are known to seek refuge in areas with structured habitat 
such as seagrass (Hovel et al. 2002; Pihl et al. 2006), the positive relationship between the 
abundance of kelp bass and different metrics (shoot frequency, shoot density and canopy 
height; Fig. 2) used to evaluate Zostera structure was expected. In California, eelgrass mit- 
igation policies established by NOAA Fisheries require that measurements of leaf shoot 
density and areal coverage (frequency, in this study) are used as a proxy for estimating 
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Table 2. Monthly estimates of total length (TL) and secondary production recorded for juvenile kelp 
bass P. clathratus in Big Fisherman’s Cove. TL values are means + SE. Zero TL values indicate survey dates 
in which individual kelp bass were either absent from the Zostera study area or >9.17 cm TL (subadult 
stage and older). 


Survey date TL (cm) =* mow!) 
Sept 2013 7.78 + 0.39 0.26 
Oct 2013 4.89 + 0.68 0.11 
Nov 2013 3.97 + 0.24 0.1 
Dec 2013 4.54 + 0.15 0.22 
Jan 2014 6.28 + 2.40 0.04 
Feb 2014 0 0 
Mar 2014 Se) se 0.37/ 0.14 
Apr 2014 0 0 
May 2014 12.0 + 0.33 0.02 
June 2014 07/3) se OZ 0.04 
July 2014 V27) a2 V533 0.45 
Aug 2014 5.87 + 0.45 0.22 
Sept 2014 5.04 + 0.09 0.57 
Oct 2014 SI) se OMS 0.44 
Nov 2014 4.36 + 0.07 0.34 
Dec 2014 4.32+0.12 0.13 
Jan 2015 0 0 
Feb 2015 5,2 se O20 0 
Mar 2015 7.50 + 0.01 0 
Apr 2015 4.33 + 0.48 0.03 
May 2015 0 0 
June 2015 9.39 + 0.43 0.13 
July 2015 BMD se OLS 0.59 
Aug 2015 OND se i (0),333) 


Secondary production 


the function of eelgrass habitat. However, a significant finding from the current study is 
that the canopy height of leaf shoots, by itself, is a better predictor of the use of Zostera 
habitat by kelp bass (Table 1). Similar sets of measurements on the spatial patterns and 
variability of seagrass-fish assemblages report that canopy height is a key determinant of 
juvenile fish density (Gullstrom et al. 2008), and leaf shoot density, on its own, is a poor 
indicator of the abundance of fishes and decapods associated with a given eelgrass bed 
(Worthington et al. 1992). In addition to providing essential habitat and protection to fish- 
ery species, seagrasses provide a variety of ecosystem services whose economic value, in 
terms of the abundance and density of seagrass beds as a whole, far outweigh their eco- 
logical function (Barbier et al. 2011; Schubert et al. 2015). Survey metrics that capture the 
three-dimensional complexity of Zostera habitat will help to better understand the ecolog- 
ical functions of these communities. 

Higher rates of secondary production in juvenile kelp bass coincided with increases in 
Zostera habitat structure (Fig. 4), as does the spawning season for this fish species, which 
reaches its peak during the summer months from June to August (Love et al. 1996; Cordes 
and Allen 1997; Erisman and Allen 2006). These findings highlight the need for time series 
data with seasonal patterns across multiple years to provide a more accurate estimate of 
variability. A reduction in the sampling frequency of fish at a given location, as dictated 
by limited access or resource availability, could yield dramatically different estimates of 
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secondary production on an annual scale. For example, if monthly rates of secondary pro- 
duction for juvenile kelp bass (Table 2) were calculated on a quarterly basis, rates would 
be either over- or underestimated by as much as 8 to 27%, respectively (quarterly sched- 
ule begins on either September or November 2013; average = 2.24 + 0.97 and 1.52 + 
0.96 g m~ yr_|, respectively; compare to 24-month average = 2.08 + 0.46 g m~ yr7!). 
Future studies should utilize both existing patterns of habitat use and spatial variability as 
a best estimate of the abundance and distribution of animals in a given location. In the case 
that fish movement patterns are limited to a specific area (cf. blue rockfish Sebastes mysti- 
nus: Jorgensen et al. 2006), a single sampling event on an annual basis may be sufficient. 

Annual rates of secondary production calculated for juvenile kelp bass in this study 
approximate those reported for adult fishes in seagrass habitat (0.24 to 7.08 g m~ yr7!, 
Table 2; compare to 3.8 to 14.6 g m7 yr_!: Adams 1976; Edgar and Shaw, 1995; Jeong 
et al. 2009). Given the increased protection of kelp bass off Catalina Island through en- 
hanced conservation and management plans (essential fish habitat and marine protected 
areas), coupled with their relatively small home range (< 0.33 ha, Lowe et al. 2003), en- 
hanced settlement and recruitment within coastal Zostera beds are likely to play an impor- 
tant role in the export of adult fish to adjacent areas (Domeier 2004; Watson et al. 2010). 
Claisse et al. (2014) report markedly higher rates of annual secondary production amongst 
fishes on offshore oil platforms (range = 104.7 to 886.8 g m~? yr!) in the same region 
as this study. Although both habitats may perform a nursery function (Love et al. 2006), 
the complex, hardscape structure of the offshore oil platform itself is distributed through- 
out the water column and 1s therefore better suited for intercepting larvae. Zostera beds, 
on the other hand, while associated with low-relief (<1 m above the benthos) nearshore 
habitats, are an important nursery ground and refuge for a variety of animals and may 
undergo high rates of species turnover for both reef fishes and invertebrates (Seitz and Ew- 
ers Lewis 2018). Nonetheless, research findings presented here indicate that Zostera is an 
important contributor to secondary production in kelp bass, which is one of the most im- 
portant nearshore recreational species off the Southern California coast (Love et al. 1996; 
Erisman and Allen 2006; Horning 2009). 

A common goal of resource monitoring is the development of biological indicators that 
provide rapid assessment of ecosystem function. Fish are useful as bioindicators because 
they are relatively easy to see, can disperse from stressed habitats and are economically 
important (Whitfield and Elliot 2002). The resources required to monitor secondary pro- 
duction at a level appropriate to detect significant impacts is daunting, but may be com- 
parable, or perhaps less intensive, than other bioindicators (cf. carbon-to-nitrogen ratios, 
rhizome sugars and above-ground biomass; McMahon et al. 2013) used to assess the qual- 
ity of Zostera habitats and how they might change over time. Further research is required 
to assess changes in secondary production in response to different impact types such as, 
sediment disturbance and invasive species encroachment. Kelp bass make for an effective 
case study of ecosystem function and services because these fish are common throughout 
the Southern California Bight and are an important component of biomass and recruit- 
ment associated with seasonal changes in Zostera habitat structure. Results from this study 
provide a baseline for the use of fish secondary production as an indicator of the health of 
eelgrass ecosystems. 

The seasonal use of coastal Zostera habitat as a nursery ground by juvenile kelp bass off 
Catalina Island during the summer and fall months coupled with the presence of signifi- 
cantly fewer fish in the winter and spring (Fig. 3 and 4) suggests that a substantial amount 
of biomass could be exported to nearby coastal habitats. These findings support previous 
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work on the export of energy and nutrients across coastal habitat boundaries through the 
transfer of fish biomass (Deegan 1993). We postulate this on the assumption that juvenile 
kelp bass vacate eelgrass beds after reaching a certain ontogenetic stage, but further work 
is needed to verify the movement of kelp bass throughout their life cycle. While the rela- 
tionship among juvenile habitats is well defined in tropical areas (Nagelkerken et al. 2000), 
far less is known in temperate regions such as the Pacific Coast. This study is one of the 
first in the region to provide quantitative data that support a possible ontogenetic shift in 
juvenile fish from coastal Zostera beds to adjacent habitats. Future studies on the export 
of fish biomass from seagrass beds (Gillanders 2006) will allow conservation and manage- 
ment efforts to focus on key marine resources such as protected areas around the California 
Channel Islands, which provide a major contribution to adult populations of ecologically 
and economically valuable species. Likewise, quantitative assessments of habitat function 
are beneficial for resource managers to determine the mitigation costs (e.g., habitat equiv- 
alency analysis) associated with a loss in natural resource services (Dunford et al. 2004). 
Further investigations of fish secondary production in nearshore coastal ecosystems are 
required to provide new insights into the mechanisms that contribute to the growth, dis- 
persal and connectivity among populations and will help to communicate the significance 
of marine resource conservation with both the general public and stakeholder groups. 
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Abstract.—Record drought from 2012 to 2016 followed by rainfall in the winter of 
2017 provided an opportunity to examine how changing climate conditions may affect 
migration opportunities for the endangered southern steelhead trout (Oncorhynchus 
mykiss). This study examined how intermittently open estuary-ocean interfaces in the 
Santa Monica Bay that have historically supported steelhead evolved temporally and 
volumetrically. All seven lagoons in the study area breached by January 2017 after 
five years of drought and nearly exclusively closed conditions. Duration of breach was 
affected by the size of the lagoon, with smaller lagoons remaining breached longer 
than larger lagoons. Conversely, volume capacity persisted longer in larger lagoons. 
Lagoon condition was quantified by presence/absence of breach and passibility, cou- 
pled with daily rainfall. This study provides important lagoon planning, restoration 
and management information needed to support recovery of southern steelhead trout 
populations in the face of climate change. 


Southern steelhead trout, Oncorhynchus mykiss, are a Distinct Population Segment 
(DPS) of west coast steelhead listed under the Federal Endangered Species Act (ESA) 
in August 1997. Historic stock numbers have drastically declined due to factors such as 
loss of estuarine habitat, land use practices that impact watershed function, and connec- 
tivity. More recently, the effect of extreme drought! that has affected rainfall and freshwater 
ecosystems across California from 2012 to 2016 has caused even more stress on these pop- 
ulations (NMFS 2016). Fewer than ten anadromous adults have been observed annually 
in the Southern California Distinct Population Segment region extending from San Luis 
Obispo to the Mexican border from 2014-2018. 

Restoration projects across the range of historic steelhead habitat in southern Califor- 
nia, and in the Santa Monica Bay in particular, have helped remove documented passage 
barriers for migrating steelhead in recent years. It should be noted that passability as de- 
fined in this study refers specifically to the lagoon — ocean connection, whereas the ability 
of fish to move upstream from the lagoon through the mainstems of creeks is usually of 
much shorter duration and more difficult to document. Three passage barriers, including a 
check dam and two Arizona crossings, were removed over the course of three years in Ar- 
royo Sequit Creek in Leo Carrillo State Park. This restoration was completed just before 


* Corresponding author: rosidagit@gmail.com 
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the rains of the 2016/2017 water year and resulted in the immigration of two anadromous 
steelhead trout (350 mm and 360 mm FL), the first anadromous individuals documented 
to return to the upper reaches of the creek since surveying began in 2004 (RCDSMM un- 
published data). Similarly, the restoration of Lower Topanga Creek in 2009 and Malibu 
Lagoon in 2012-2013 have resulted in increased critical fish passage opportunities allowing 
three anadromous fish to migrate upstream from the ocean with the 2017 rains (RCDSMM 
unpublished data). 

Despite the success of restoration and anthropogenic passage barrier removal projects, 
natural passage barriers still exist for O. mykiss due to the nature of intermittently open 
bar-built coastal lagoons that are characteristic of Mediterranean climate regions, such as 
those in the Santa Monica Bay (McLaughlin et al. 2014). These “small creek” archetypal 
systems are identified as intermittently open bar-built estuaries.* According to the lagoon 
classification system by Jacobs et al. (2011)*, the small, steep coastal creeks draining into 
the Santa Monica Bay respond similarly to the complex interactions of wave exposure, 
longshore transport, sediment delivery, watershed size, and formation processes. While his- 
toric data on creek mouths provided in the US Coast and Geodetic Topographic Survey 
(T-sheets) provides interesting information on past conditions, we found that due to exten- 
sive modifications over time, bar-built estuaries along the Santa Monica Bay needed to be 
further examined to better understand present constraints. 

Lagoons not only provide important fish passage opportunities, but also can provide 
important foraging habitat for smolts - young trout migrating outward to the ocean (Bond 
et al. 2008). In central and northern California systems, the opportunity to spend time in 
brackish water adjusting to higher salinity levels also lowers stress levels for smolts dur- 
ing the freshwater to ocean transition (MacDonald et al. 1988). In order to support the 
recovery of southern steelhead trout, adequate connectivity is needed for immigration of 
anadromous adults to spawn and for juvenile smolts to migrate out to the ocean. If this 
exchange of habitat is unable to occur, gene flow and fecundity will continue to be reduced 
in the populations within the Santa Monica Bay, making them more susceptible to extir- 
pation (NMFS 2012). The combination of current anthropogenic constraints caused by 
development adjacent to the creek mouths, impacts of Pacific Coast Highway and beach 
bar manipulation, muddle the behavior one would expect from examining past data and 
color restoration opportunities and constraints. This information will help us to under- 
stand how bar built estuaries in the Santa Monica Bay might respond to future climate 
change and sea level rise. 

Systematic studies of the processes driving estuarine seasonal closure dynamics through- 
out Mediterranean zones worldwide provide some guidance, although observations of 
the patterns in southern California appear to be somewhat different than those observed 
elsewhere. McSweeney et al. (2017) developed a classification system for intermittently 
open/closed estuaries (IOCE) in Victoria, Australia, a Mediterranean climate like southern 
California. However, lagoons in Australia are tidally influenced, much larger in all ways, 
and unconstrained by topography so their behavior patterns differ. Northern California 


* Southern California Wetlands Recovery Project. 2018. Wetlands on the edge: the future of south- 
ern California’s wetlands: Regional Strategy 2018. Prepared by the California State Coastal Conservancy, 
Oakland, CA. 

* Jacobs, D., E. D. Stein, and T. Longcore. 2011. Classification of California estuaries based on natural 
closure patterns: Templates for restoration and management. Southern California Coastal Water Research 
Project Tech. Memo. 619:1-50. 
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lagoons are much more exposed to wind and wave action due to their west facing orienta- 
tion, yet many of the lagoons examined in northern California also breach naturally under 
lagoon-side flow driven processes instead of tidal processes (Kraus et al. 2008). A South 
African definition for estuaries also echoes processes seen in our study area. Describing a 
lagoon as “a coastal body of water in intermittent contact with the open sea and within 
which sea water is measurably diluted with fresh water from land drainage” (Day 1980), 
takes into account systems that may be dry for periods of time (Cooper 2001). This trait 
is seen in many of the lagoons in the Santa Monica Bay that dry down completely during 
summer/fall periods when rainfall is typically absent. Even when lagoon surface area per- 
sists through the dry season, connection between a lagoon and upstream habitat is often 
cut off beginning in the early summer due to subsurface flow in lower reaches (Dagit et al. 
2017). 

Due to the low coastal exposure of the south/southwest facing beaches, lagoons of the 
Santa Monica Bay are not strongly influenced by wave action.t The Santa Monica Bay 
coastal area is situated facing south/southwest, which reduces the importance of wind 
action on beaches by diminishing strong wave action that contributes to longshore sed- 
iment processes that help deteriorate berms in higher wave exposure areas.” Dominant 
wave action is due to prevailing northwest winds during spring and summer months, with 
westerly winds that are common in fall through spring. However, south-southeast winter 
storm winds can sometimes generate strong waves that cause significant erosion. Overall, 
wave climate in the area is considered mild due to sheltering from deep ocean waves by the 
California Channel Islands, although west-facing beaches are more exposed to swell from 
storms in the North Pacific, and south facing beaches more susceptible to waves generated 
by southern swell. The lagoons east of Point Dume (including Solstice, Las Flores and 
Topanga), are generally more sheltered than those west of Point Dume (Zuma, Trancas, 
Arroyo Sequit and Big Sycamore). Waves greater than 3.5 m are most likely to erode sand 
from the beach, but these waves have not been found to occur often.® Instead, stream in- 
put and rainfall are the determining factor of breach potential and surface water presence. 
Stream flow, which is highly seasonal in California, is a major control on the state of a 
lagoon’s inlet (Elwany et al. 1998).* Rich and Keller (2012) examined how watershed size 
and upstream slope can predict the behavior of small lagoons along the Santa Barbara 
Coast, determining that watershed size and the size of the associated lagoon are positively 
correlated. 

Additionally, all the lagoons in the study zone are considered hydraulic estuaries that are 
normally closed to the sea during drier months by a bar that forms across the mouth of the 
lagoon.* Higher storm generated waves and higher tides both occur during the winter and, 
combined with higher flows associated with storm events, are associated with breaching 
the sand berms and re-connecting coastal creeks and lagoons with the ocean. Another 
characteristic that we did not specifically examine is tidal prism, the volume of water in a 
lagoon that is between mean high tide and mean low tide. Tidal prisms in river-dominated 
lagoons, where estuary side processes initiate breaches, are too small to maintain a breach 
against the wave action and tidal influence that the beach faces in the Santa Monica Bay 


* Moffatt and Nichol. 2013. Broad Beach Restoration Project Coastal Engineering Report Exhibit L to 
CDP Application 4-12-043. Prepared for Broad Beach Geologic Hazard Abatement District. 

° Noble Consultants. 2009. Coast of California storm and tidal waves study, Los Angeles Region. Pre- 
pared for U. S. Army Corps of Engineers. 
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typically receive (Cooper 2001). Like coastal exposure and formation process, tidal prism 
characteristics in the study area lagoons are defined more by morphological constraints 
that are uniform throughout the study area; additionally tidal prism is usually so small 
that it is not able to initiate breaches (Cooper 2001). 

This study examined how intermittently open bar-built lagoon morphology, specifically 
breach frequency and volumetric changes, responded to rainfall after a period of drought, 
and the implications this may have on future restoration efforts to improve O. mykiss pas- 
sage opportunities in the Santa Monica Bay. Seven lagoons within the Santa Monica Bay 
were selected for evaluation. Historically, Arroyo Sequit, Big Sycamore, Las Flores, Mal- 
ibu, Solstice, Topanga, Trancas, and Zuma Creeks all supported southern steelhead trout 
populations (NMFS 2012), but current lagoon conditions restrict passage opportunities at 
all but Malibu Creek, which was restored in 2012. Malibu Creek is the only antecedent 
creek that crosses the Santa Monica Mountains from the San Fernando Valley to the 
ocean; its course remained the same and cut through bedrock despite the uprising of the 
Santa Monica Mountains. Malibu Creek also receives extensive imported water flows. Due 
to the recent restoration actions undertaken at this lagoon and its differences from the 
other coastal creeks, which do not experience either significant inputs of imported water 
or withdrawal of groundwater resources, Malibu is not further considered in this study. 
The baseline in-depth data provided for the seven lagoons in this study provide critical 
context for land use managers and restoration planners as the systems in our study have 
been identified as a priority watersheds for restoration (NMFS 2012). 


Materials and Methods 


Based on the watersheds identified in the Southern California Steelhead Recovery Plan 
(NMES 2012), the following lagoons and creeks within the Santa Monica Bay were the 
focus of this study: Arroyo Sequit, Big Sycamore, Las Flores, Solstice, Topanga, Trancas, 
and Zuma. Although Malibu Creek is the second largest watershed draining into the Santa 
Monica Bay with present and historic southern steelhead trout habitat, Malibu lagoon was 
not included in this analysis. Due to the augmentation of flow from Tapia Water Reclama- 
tion Facility’s discharge in accordance with California Regional Water Quality Control 
Board Los Angeles Region’s NPDES NO. CA0056014, Malibu Creek and lagoon are the 
only system within the Santa Monica Bay that has a significantly altered hydrologic profile. 
All the other lagoons and their associated creeks rely on streamflow provided by natural 
variables such as rainfali and/or the water table. Regional context of these watersheds 
within the Santa Monica Mountains National Recreation Area and Santa Monica Bay 
(Fig. 1), along with individual watershed boundaries and limits of anadromy (Fig. 2) for all 
the study watersheds and lagoons are important to note when considering anthropogenic 
alterations. All of the systems in our study area have been altered in some way, chang- 
ing the length of stream available for southern steelhead trout and imposing lagoon con- 
straints due to Pacific Coast Highway (PCH) bridges (Table 1). Other anthropogenic fac- 
tors include channelization in creek reaches, concrete culverts, multiple bridge constraints, 
Arizona type instream crossings, and fully restricted lagoons by wall or riprap barriers. 

Rainfall metrics were obtained using LA County Department of Public Works, Water 
Resources Division’s “Precipitation” webpage (http://www.ladpw.org/wrd/precip/). Us- 
ing the “Near Real-Time Precipitation Map”, and gauges with closest proximity, daily 
rainfall was calculated for each lagoon’s watershed. The Lechuza gauge (# 454) data was 
used for Arroyo Sequit, Big Sycamore, Trancas and Zuma watersheds. The Big Rock Mesa 
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Service Layer Credits: Sources: Esn, USGS, NOAA 


Fig. 1. Watershed boundaries of the northern Santa Monica Bay. 


gauge (#320) was used for Las Flores watershed, the Monte Nido gauge (#319) was used 
for Solstice watershed, and the Topanga Canyon gauge (#318) was used for Topanga wa- 
tershed. Flow data was obtained for Topanga Creek from the Topanga surface water flow 
gauge (F54C) maintained by Los Angeles County. Data was available for 2013-2017, but 
due to changes in thalweg location, the gauge currently sits a couple meters away from 
the thalweg and does not capture any low flow measurements, despite flow being present. 
Although we recognize the importance of flow, because this data was so limited and not 
available for all creeks, it was not incorporated further in this study. 

Overwash from high tides and storm event tides can affect passage opportunities for 
southern steelhead trout into the study area’s lagoons. Tide data for the 2013-2017 study 
period was measured at Los Angeles Station ID: 9410660, (33.7200, -118.2720). Because 


Table 1. Descriptive measurements of each lagoon and watershed. 


PCH Hwy Mean 

Watershed/ Limit of Culvert Annual Coastal Average Lagoon 

Location Catchment Area Anadromy Dimensions Rainfall’ Exposure Wetted SA, 2017 
Arroyo Sequit 31km? 2.64rkm' 24.2mx23.75m 14.094 Low NOB wie 
Big Sycamore 54.4km? 19.3rkm 16.8m x 23.2m 14.094 Low 2009.67 m? 
Las Flores 11.7km? 3.2rkm 16.6m x 27.7m 11.874 Low 557 m? 
Solstice 11.1km? 30a  PjNiax4syan NSA Low 140.33 m? 
Topanga 46.6km? 5.3rkm 19.7m x 27.5m i Low 4561.33 m? 
Trancas 26.4km? 10.46rkm  =30.5m x 25.8m 14.094 Low 1995.67 m? 
Zuma 25.9km? 8.05rkm  30.4m x 23.4m 14.094 Low 9624.33 m? 


* Mean annual rainfall calculated by averaging rainfall from 2013-2017 water year data collected. 
‘Limit of anadromy under low flow conditions. 
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Fig. 2. Individual Watershed boundaries and creek locations. 


we describe breach frequency when the connection was both passable and not passable, 
passable conditions were augmented for times when high tide could allow passage. When 
a lagoon was breached, any tides over 1.52 m above Mean Low Low Water (MLLW), were 
noted for potentially passable connection during high tide. 

~ Four comprehensive longitudinal profile lagoon surveys occurred in November 2013 
and January, June, and October 2017 at all sites. Lengths and widths of wetted area were 
recorded at wetted cross sections using hand held 100-meter tapes. Garmin GPS units 
recorded waypoints of lagoon boundaries (DATUM: WGS84) and GPS locations were up- 
loaded into Google Earth. Lagoon area polygons were then mapped using the waypoints, 
observable landmarks (such as the PCH bridges), and survey notes. Storm-event related 
monitoring of the lagoon/ocean interface was conducted from October 2013 through De- 
cember 2016 in order to document passage opportunities and constraints. In response 
to large and frequent storms starting in January 2017, weekly lagoon monitoring was 
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conducted by a field team of two or three surveyors. These surveys consisted of visual doc- 
umentation of lagoon and breach conditions via field notes and photo points to record 
physical condition (depth, connectivity, etc.) changes over time. Passibility was also 
recorded at the time of these surveys. Any connection between the lagoon and the ocean, 
even a steady connected trickle, was considered a breach. Only when flow was of adequate 
depth and size was the breach considered and noted as passable. A connection was consid- 
ered passable if there was at least 15-30 cm of depth in the connection between the ocean 
and lagoon.’ It is possible for a connection to be passable only at high tide, which was dif- 
ficult to assess 1n our Surveys since each site was assessed at varying times throughout the 
tidal cycle. Subsequently, a lagoon was only marked passable if it was clearly accessible at 
the time of survey. Later tide analysis was performed to predict additional potentially pass- 
able opportunities at high tide, a roughly 1-3 hr window. When a lagoon was breached, any 
tides over 1.52 m (MLLW), were examined for passable connection during high tide. Data 
from less frequent breach condition surveys from 2013-2016 was included for temporal 
context. 

Field measurements of the lagoons documented conditions in January, June, and 
October 2017; at the beginning of rains, the beginning of summer, and the end of the water 
year when the lagoons have been subjected to the entire dry season. Physical measure- 
ments were then mapped in Google Earth to create a visual time series of lagoon dimen- 
sions and to obtain temporal descriptive measurements of lagoon wetted surface area (SA) 
and volume (V). SA and V from November 2013 was included in tables for context but 
not included in analysis since survey methods and measurements were not consistent with 
the surveys conducted in 2017. Average depths were estimated based on multiple depth 
measurements taken with a meter stick throughout each lagoon, but this method could 
be improved by using water level data loggers that would give a continuous timeline of 
changing lagoon depth over time. This method was not utilized in this study due to limited 
funding for such loggers. 

All data was entered by date and location. Graphs of lagoon breach condition in relation 
to rainfall on a temporal scale were created for each lagoon. Lagoon volume was calculated 
by mapping lagoon boundaries via Google Earth mapping software, creating polygons to 
the specific measurements calculated and recorded in the field. Surface area calculations 
were coupled with average depth measurements recorded during survey visits to calculate 
an estimated volume at that particular point in time. An observed rate of rainfall required 
to cause a breach was calculated by taking the total amount of rain fallen in the watershed 
from October 1*' to first day of breach for each lagoon and dividing it by the number of 
days until first breach, SS ae Statistical analysis was performed using JMP 
Statistical software. Continued surveys following the same methods may provide a data set 
large enough to perform robust statistical analysis in the future. 


Results 


In the 2016-2017 water year, rainfall throughout the Santa Monica Bay averaged 
65.79 cm. Rain data for each watershed can be found in Table 2. Rain data was used 
to calculate several variables that were used for volume and breach frequency analysis. 


Flosi, G., and F. L. Reynolds. 2010. California salmonid stream habitat restoration manual. Second edi- 
tion. California Department of Fish and Game, Sacramento, CA. http://www.dfg.ca.gov/fish/Resources/ 
HabitatManual.asp. 


180 SOUTHERN CALIFORNIA ACADEMY OF SCIENCES 


Table 2. Watershed and rain characteristics of study area for the 2017 water year. 


Watershed Area 2017 Rainfall Watershed Rain Rain Rate 

Location (km?) (cm) Volume (m+?) (cm/day) 
Arroyo Sequit 31 73.81 22,881,800 0.092 
Big Sycamore 54.4 73.81 40,153,900 0.162 
Las Flores il 56.21 6,576,590 0.076 
Solstice lI 62.9 6,992,110 0.069 
Topanga 46.6 67.03 31,236,300 0.056 
Trancas 26.4 73.81 19,486,500 0.082 


Zuma BJS) 73.81 19,117,400 0.148 


These include the volume of rain fallen (watershed area multiplied by rainfall) over each 
watershed during the 2017 water year and rain rate, defined in this study as centimeters of 
rain per day until first breach. 

Topanga Creek is the only system that has a flow gauge, but it is currently located in a 
position that does not accurately measure low flow conditions. According to the available 
data, Topanga Creek had a mean daily discharge of 0.11 m/s, a maximum of 1.66 m°?/s, 
and a minimum of 0 m?/s for the 2017 water year. Average daily discharge was 0.024 m?/s 
for the 2012 water year, 0.021 m*/s for 2013, 0.017 m?/s for 2014, 0.027 m?/s for 2015, and 
0.023 m?/s for 2016. No flow data was available for any other site. A conceptual hydrology 
analysis estimated flow rate in cubic feet per second for Trancas lagoon for several po- 
tential lagoon restoration alternative designs. Depending on lagoon restoration area size, 
determined that a two-year storm could produce flows between 0.007-0.033 m3/s, while a 
ten year storm could produce flows between 0.028-0.282 m?/s.° 

All lagoons showed a time lag between end of breach and end of observed overwash, an 
expected observation as sediment deposition to rebuild a berm takes several days to weeks, 
depending on the specific location. Seaward wave action, already of low influence for the 
study area, became less of a factor on breach potential when a berm barrier was high and 
wide. When accounting for tidal influence, the amount of time of passable connection was 
increased by tides when a system was breached. Recorded days of full passibility for each 
lagoon can be found in Table 4. Due to funding constraints, it was not possible to measure 
tidal prism as part of this study. 

All the lagoons — with the exception of Topanga — were consistently dried down almost 
completely between 2013 and 2016, highlighting the severity of the California drought 
(Table 3). Arroyo Sequit had a modest volume throughout the 2017 water year with a max- 
imum of 821 m? in June 2017. Volume in Big Sycamore measured less than 100 m? during 
both November 2013 and January 2017 survey, and then increased more than a factor of 
20 by June 2017. Anthropogenic constraints limited any lateral movement of water in Las 
Flores and resulted in a maximum volume of only 337 m°. Similarly, constraints on Solstice 
Lagoon from the PCH culvert and manipulation by individuals (placing sand bag, digging 
small trenches) resulted in a small maximum volume of 125 m3. Trancas Lagoon had a 
dramatic volumetric change, quadrupling in wetted surface area and doubling in estimated 
depth throughout the study period. Topanga lagoon held onto its surface water the longest 


i Dagit, R., S. Albers, and C. Stevens. 2015. Trancas lagoon restoration feasibility study 2013-2015. Pre- 
pared for CDFW and Zuma Beach Properties, LLC. Resource Conservation District of the Santa Monica 
Mountains, Topanga. 
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Table 3. Recorded lagoon surface area, depth, and volume including changes in measurements between 


November 2013, January 2017, June 2017, and October 2017 surveys. 


Location 


Arroyo Sequit 


Big Sycamore 


Las Flores 


Solstice 


Trancas 


Topanga 


Zuma 


Date of 
Survey 


11/1/2013 
11/1/2014 
11/1/2015 
11/1/2016 
1/27/2017 
6/23/2017 
10/1/2017 
11/1/2013 
11/1/2014 
11/17/2015 
11/1/2016 
1/27/2017 
6/23/2017 
10/1/2017 
11/1/2013 
11/1/2014 
11/1/2015 
11/1/2016 
1/27/2017 
6/23/2017 
10/1/2017 
11/1/2013 
11/1/2014 
11/1/2015 
11/1/2016 
1/27/2017 
6/23/2017 
10/1/2017 
11/1/2013 
11/1/2014 
11/1/2015 
11/1/2016 
1/27/2017 
6/23/2017 
10/1/2017 
11/1/2013 
11/1/2014 
11/1/2015 
11/1/2016 
1/27/2017 
6/23/2017 
10/1/2017 
11/1/2013 
11/1/2014 
11/1/2015 
11/1/2016 
WATE) 
6/23/2017 
10/1/2017 


Wetted 
Area 


A Wetted 
Area 


No Data 
No Data 
No Data 
No Data 
1036 
-1282 


-813 


Average 


Depth 


No Data 
No Data 
No Data 


0.1 

0.75 
0.75 
0.16 


S 
S 
rs 


eS) 
i) 


SSSeeSeeeees 
oS) 


wo — 
onion 


No Data 

No Data 

No Data 
-0.65 
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into the drought with some surface water consistently present (no surface area/volume 
measurements available) through 2014-2016 and an increase in volume measuring 4000 m? 
during the June 2017 survey. Zuma lagoon was the quickest to recover from drought condi- 
tions, drying down almost completely from 2014-2016 and then increasing to an estimated 
volume of 15,222 m? in the January 2017 survey. 

Throughout the 2017 water year, a trend in volume changes was observed. Apart from 
Big Sycamore, all lagoons had a small volume in January 2017 when all lagoons were 
breached and outflowing, recorded their largest volume measurement in June 2017, and 
measured the smallest volume in October 2017 because of dry down from high summer 
temperatures. Big Sycamore was the only location that had a larger volume in October 2017 
than in January 2017. Lagoon depth was extremely variable throughout the study area. 
Some locations increased in depth between each survey in 2017 (Trancas and Topanga); 
some locations decreased in depth between each survey (Arroyo Sequit and Zuma); and 
some locations increased depth between January and June, then decreased in depth be- 
tween June to October (Big Sycamore, Las Flores, Solstice). Fig. 3 and Fig. 4 illustrate the 
changes in lagoon volume over time at each location in this study. Some mapped surface 
areas in Fig. 3 appear to extend into the surf. GPS data for these instances were verified 
for accuracy. Though all surveys were conducted during the day in the same order, some 
occurred during low tide times, which explains the appearance of some of the maps ex- 
tending into the ocean. 27 January 2017 was the day of the new moon so neap tides were 
experienced with an afternoon low tide of 0.027m below MLLW. 

Anthropogenic factors affect all the study locations and influence the volume capacity 
in some way. Instream creek barriers in Arroyo Sequit were recently restored, but the PCH 
bridge, a concrete pedestrian walkway, and a vehicle driveway line the edges of the lagoon, 
confining its footprint. Big Sycamore is under similar constraints with a pedestrian bridge 
running across the upper portion of the lagoon and the PCH bridge along the lower end 
where riprap holds the banks in place between these two bridges. Las Flores is bordered 
on either side by tall concrete retention walls with buildings on the east bank, and Solstice 
has a long culvert under PCH which causes sheet flow into the lagoon and cuts off any 
upstream migration possibility. Topanga also suffers from constraints due to the concrete 
retention walls extending from the PCH bridge and fill slopes to a regularly groomed berm. 
Unlike the other lagoons with more natural channels upstream of the mouth, Trancas 
Creek has a 600 m concrete trapezoidal flood control channel that begins upstream of the 
lagoon and channels flow in a way that may affect lagoon morphology. Zuma lagoon has 
the most unconstrained banks and is by far the largest in our study, but multiple Arizona 
crossings, regular grooming to maintain the berm for lifeguard vehicle access, and vector 
control manipulation within the lagoon itself all affect its suitability as potential O. mykiss 
habitat. 

Increasing lagoon volume correlated with higher surface areas in the study area lagoons 
. (R? = 0.897, p < 0.001). Due to the multi-year drought, all lagoons effectively started the 
2016-2017 water year almost completely dry, but the ensuing 2017 rainy season increased 
surface flows and restored passable conditions into the upper reaches of each system. Santa 
Monica Bay lagoons with bigger surface areas were more likely to contain a larger volume 
of water (Fig. 5). There was no significant relationship between watershed area and lagoon 
surface area. Controlling for differences in rain received over different watersheds in the 
study area, total rainfall for the 2017 water year was couple with watershed area to get 
a volume of rain fallen within each watershed (Table 2). Watershed volume in relation to 
lagoon surface area yielded no significant relationship. A more robust sample size in the 
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Fig. 3. Mapped surface area of Arroyo Sequit, Big Sycamore, and Las Flores lagoons on surveys con- 
ducted on 3 November 2013, 27 January 2017, 23 June 2017, and 5 October 2017 (maps ordered chrono- 
logically from left to right). *Las Flores Lagoon dried down completely when observed a week prior to the 
5 Octobe4 2017 survey. County employees bulldozed the lagoon, removing all in-channel vegetation as well 
as several inches of sand revealing sub-surface water that was not previously exposed. 


future would improve the confidence in our analysis and be a better indicator of the true 
relationship between catchment area and lagoon size. Fig. 6 shows how depth in the study 
area lagoons evolved during a year of average rainfall following a multiyear drought. Depth 
peaked in the middle of the year, when rains had ceased but before the peak in average high 
temperatures that Southern California experiences from July to September each year.’ 

A summary of storm events, breaches, and possible fish passage opportunities is pro- 
vided in Table 4. Passable connection was defined as a breach with enough depth for 


” US Climate Data 2018, Climate — Los Angeles. https://www.usclimatedata.com/climate/los-angeles/ 
california/united-states/usca 1339. 
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Fig. 4. Mapped surface area of Solstice, Topanga, Trancas, and Zuma lagoons on surveys conducted 


on 3 November 2013, 27 January 2017, 23 June 2017, and 5 October 2017 (maps ordered chronologically 
from left to right). 
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Table 4. Summary of rainfall and fish passage opportunities of monitored lagoons in the Santa Monica 
Mountains, in alphabetical order. 


Estimated 
Rainfall Number Number of 
Total” of Days Days Passable —O. mykiss 
Lagoon Water Year _—(inches) Dates Entrance Open Breached to O. mykiss present 
Arroyo Sequit 2013-2014 5.87 unknown 0 0 Yes 
2014-2015 WES) none 0 0 Yes 
2015-2016 9.09 29 Oct 2015 l l Yes 
2016-2017 29.06 30 Dec 2016- 14 Jan 97 5] Yes 
2017 
18 Jan — 8 Apr 2017 
Big Sycamore 2013-2014 5.87 none 0 0 No 
2014-2015 NWS 03 Mar 2015 (trickle) 0 0 No 
29 October 2015 
(storm) 
2015-2016 9.09 none 0 0 No 
2016-2017 29.06 23 — 27 Jan 2017 39 30 No 
3 Feb — 8 Mar 2017 
Las Flores 2013-2014 5.16 26 Feb — 02 Mar 2014 4 4 No 
2014-2015 WA 12 Dec 2014 — 3 Mar 10 <10 No 
2015 (pinched) 
2015-2016 7.36 6 — 12 Mar 2016 7 No 
2016-2017 DIM 24 Dec 2016 — 21 Apr 119 44 No 
2017 
Solstice 2013-2014 6.31 none 0 0 No 
2014-2015 16.06 12 — 16 Dec 2014 7 <10 No 
12 Jan 2015 
2015-2016 9.84 5 — 15 Mar 2016 1] 0 No 
2016-2017 24.8 23 Dec 2016 — 3 Jul 193 5 No 
2017 
Topanga 2013-2014 6.84 27 Feb — 07 Mar 2014 10 4 Yes 
25 March 
03 Dec and 12-16 Dec 
2014 
2014-2015 13.76 10 — 19 Jan and 3 Mar 17 <10 Yes 
2015 
03 Mar 2015 
2015-2016 8.31 6 March — 6 April By) 5 Yes 
2016 
2016-2017 26.39 21 — 26 Dec 2016 132 36 Yes 
11 Jan — 16 May 2017 
Trancas 2013-2014 Dat d/ none 0 0 No 
2014-2015 17.5 3-8 Dec 2014 5 5) No 
(overwash) 
16 Dec 2014 (pinched) 
2015-2016 9.09 6 Mar 2016 l 0 No 
2016-2017 29.06 24 Dec 2016 — 12 Jan 110 61 No 
2017 
19 Jan — 18 Apr 2017 
Zuma 2013-2014 5.87 none 0 0 No 
2014-2015 IVES none 0 0 No 
2015-2016 9.09 none 0 0 No 
2016-2017 29.06 20 Jan — 1 Mar 2017 4] 39 No 


* Water year, as used by the USGS: October 1 through September 30 of the following year. 
** Rainfall data from gauge stations managed by Los Angeles County Department of Public Works. 


186 SOUTHERN CALIFORNIA ACADEMY OF SCIENCES 


Wetted Surface Area vs Estimated Volume 
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Fig. 5. Wetted surface area measurements are positively associated (p < 0.001) with increased volumes 
for all lagoon surveys conducted post-rain in 2017. 


O. mykiss to migrate to or from the ocean to the lagoon during all tides. Breach patterns 
for lagoons with larger surface areas such as Big Sycamore and Topanga were strongly 
influenced both by rain events and by the constraints of the PCH bridges and fill slopes 
that border these lagoons. Breaches closed rapidly after storms in the beginning of the 
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Fig. 6. Average depth measurements of each lagoon during 27 January 2017, 23 June 2017, and 5 Oc- 
tober 2017 survey. Includes mean depth (line curve; cm) and confidence of fit (background shading). 
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season, but once reopened, the lack of ability to spread out helped sustain subsequent 
breaches. Zuma lagoon was more influenced by rain events, as this location lacked lagoon 
bank constraints. It took less cumulative rainfall to cause a breach in Topanga and Big 
Sycamore than in Zuma, but once breached, Zuma’s connectivity to the ocean was sus- 
tained. Breaches for the lagoons with smaller footprints tended to occur more quickly dur- 
ing the rainy season and lasted for months. Small windows of closure occurred in Arroyo 
Sequit and Trancas, but they did not close for long. It should be noted that passibility in 
this study refers specifically to the lagoon—ocean connection, whereas the ability of fish to 
move upstream from the lagoon through the mainstems of the creeks is usually of much 
shorter duration and more difficult to document. Documented windows of passibility dur- 
ing a site visit, when a breach was clearly passable for O. mykiss, were recorded and later 
analysis of tides resulted in an estimation of possible passable windows during high tide to 
give a detailed evolution of breach at each location over the span of the 2017 water year 
(Fig. 7). No relationship was identified between watershed volume and number of days of 
passable connection. There was more significance to the relationship between watershed 
volume and the number of days a system was breached. 

The open period for all lagoons ranged from 39 d (10.68% of the year) to 193 d (52.88%) 
in the 2016-2017 water year, with an overall average time breached at 104 d (28.61%). 
Comparing the number of days breached with the average lagoon surface area measured 
on 27 January, 23 June, and 5 October 2017, there was a slight trend showing that larger 
surface area lagoons were open for fewer days, and smaller surface area lagoons were open 
for more days over the water year. Since the surface area values for each location were 
averaged between the three survey events, there are only seven data points to work with 
and no significance was calculated in this relationship. 

Calculated rain rates ranged from 0.143 to 0.412 cm/d with a median of 0.221 cm/d. 
Comparing the rate of rainfall/d needed for first breach by the average lagoon surface area 
as measured on 27 January, 23 June, and 5 October 2017, larger lagoons needed more rain 
to initiate first breach, whereas smaller lagoons were able to breach with less accumulated 
rainfall. When examining the effects of lagoon surface area and watershed volume on rain 
rate needed for first breach (Fig. 8), both were positively correlated. Increasing average sur- 
face area or watershed volume predicted an increase in rain rate (i.e. an increased amount 
of rain required for first breach to occur), though linear regression revealed a higher confi- 
dence in the line of best fit between rain rate and watershed volume. A significant negatively 
correlated relationship between rain rate and the number of days a system was breached 
was discovered. Locations that were breached the longest also had lower rain rates, while 
shorter breach windows required a higher rain rate (R? = 0.763 and p = 0.01). 


Discussion 


Restoration planning is in progress for many of the lagoon systems in the Santa 
Monica Bay and design for these efforts not only must account for site-specific conditions, 
but also address regional wetland restoration and species recovery goals. When considering 
restoration options, it is critical to evaluate the implications of lagoon size because breach- 
ing frequency may depend on it (Rich and Keller 2012). Data on the presence/absence 
of steelhead, as well as information regarding the existing lagoon breach patterns, was 
compiled in order to identify the appropriate regional archetypes and lagoon classification 
patterns, which can help guide restoration designs. The challenge for restoration planners is 
to provide the right balance between connectivity during the winter months when stream 
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Fig. 7. Lagoon breach condition related to rainfall over the 2017 water year. (Breach = dark grey bar; 
Passable Connection = light grey bar; Rainfall = dark line columns). 
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Fig. 7. Continued. 
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Fig. 7. Continued. 


flows typically drive lagoon breaches (McLaughlin et al. 2014)*, versus accommodating 
more wetland habitat and transition zones that will provide flexible response to projected 
and measured sea level rise (Nerem et al. 2018). When stream flow is a limiting factor, 
the general assumption that larger systems remain open for longer periods of time when 
compared to their smaller counterparts was not observed during our study.* 

All seven of the lagoon-creek systems examined in this study fit the “small creek” 
archetype, and while similar patterns were observed in regard to breach dynamics, fur- 
ther classification was confounded by the unique constraints particular to each system. 
Rich and Keller’s (2012) analysis of small coastal lagoons near Santa Barbara, CA, a sim- 
ilar coastal area, found that the size of an estuary was directly related to watershed size. 
In their study, larger watersheds produced larger creek channels which in turn resulted in 
larger lagoons. While this relationship could not be confirmed in our study area, other 
observations from locations with similar characteristics were also seen in the Santa Mon- 
ica Bay. As has been observed in several Mediterranean regions, breaching was initiated 
by processes on the lagoon side of the berm, namely the lagoon/ocean water elevation 
difference and the width of the barrier beach (Kraus et al. 2008). Berm closure was ini- 
~ tially caused by wave driven sediment import that occurred at a rate slightly higher than 
sediment export from tides and creek flow (FitzGerald 1996). Behrens et al. (2013) also 
observed that the amount of closed days is more dependent on river flow and less related 
to wave condition. These observations closely match the patterns documented in this study. 
Despite differences in volume capacity, the drought reduced the lagoons in our study area 
to meager puddles or completely dry channels that no O. mykiss were able to access for 
years. Evaporation rates and water seepage due to the porosity of substrate are factors that 
can contribute to lack of surface water in lagoons, and the lowering of the water table likely 
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Fig. 8. Rate of rain/d versus number of days breached. 


contributed to an overarching dry down that extended upstream from the lagoons during 
the drought period. 

Of significance for O. mykiss in particular, the duration of a breach often exceeded the 
window of possibility when there was sufficient depth for either returning anadromous 
adults or emigrating smolts. Adequate depth is necessary to allow for O. mykiss pas- 
sage, with a minimum 0.305 m required for anadromous adults and 0.152 m required 
for migrating juveniles (Flosi and Reynolds 2010). Sustaining adequate depth for mi- 
gration of O. mykiss both in and out of a creek system is a key factor that needs to 
be taken into consideration when restoration planning occurs for a lagoon in the Santa 
Monica Bay. 

Southern steelhead trout life history cycles require three important habitat conditions: 
connection from the ocean into the lagoon, suitable lagoon conditions for smolts to tran- 
sition from freshwater to saltwater (Bond et al. 2008), and connection from the lagoon 
upstream to upper creek habitats allowing passage for both smolts and anadromous adults 
(NMES 2012). Recovery of this species relies upon all of these elements being functional 
simultaneously in numerous coastal streams. Connectivity both to the ocean and upstream 
was extremely limited during the study. Depths measured in these lagoons were on the 
lower limit of the preferred ranges for supporting smolt growth or transition within the 
lagoon itself, and the limited wetland and transitional habitat in most of these locations 
provide limited support for smolts. 

The observed ability of larger lagoons in our study area to retain water through the sum- 
mer and into the following water year has important implications for O. mykiss smolts. 
Smaller smolts with an average fork length of 112 mm were found to stay in a central coast 
estuary and wait until the following winter to emigrate when they were much larger (Bond 
et al. 2008). Smolts would not have this opportunity in the smaller lagoons in the Santa 
Monica Bay, as Arroyo Sequit, Solstice, and Las Flores lagoons dried down to uninhabit- 
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able levels as the study period progressed into summer 2017. Data examining the average 
fork length of smolts in Topanga Creek indicate that most emigrating smolts are larger 
than 170 mm (Bell et al. 2011), which suggests that environmental cues in the Santa Monica 
Bay may favor larger smolt sizes for emigration. Even with adequate surface water remain- 
ing present throughout the summer in lagoons like Big Sycamore, Topanga, Trancas, and 
Zuma, connection to upstream habitat was cut off very early in the dry season, between 
May and June. This phenomenon had been observed consistently in Topanga Creek for 
many years (Dagit et al. 2017). The lower reach of Topanga Creek has been drying down 
earlier each successive year of drought, with the length of reach dry also increasing each 
year. 

Previous habitat assessment and barrier identification work in the Santa Monica Bay 
documented watershed size and geomorphology, both natural and anthropogenic fish pas- 
sage barriers and habitat quality (CalTrout 2006). Based on an integrated assessment 
of these factors, CalTrout (2006) prioritized restoration actions to focus on recovery of 
O. mykiss. Their recommended ranking prioritized Topanga, Arroyo Sequit, Trancas, Big 
Sycamore, Zuma, Las Flores and Solstice. Our study builds on that effort by providing 
both qualitative and quantitative data on the evolution of breaching and lagoon condi- 
tions. Based on the results of this study, the unique challenges and constraints of each sys- 
tem are discussed, with ranking based on all of the above metrics, as well as documented 
presence of O. mykiss since 2000 (RCDSMM unpublished data). 

Because it has the best average habitat quality and the third largest habitat quantity in 
kilometers, as well as supporting the only reproducing population of O. mykiss remaining 
in the Santa Monica Bay, restoration of Topanga lagoon is considered to be the highest 
restoration priority. Topanga lagoon’s unique system sustains a steady pool of surface wa- 
ter throughout the year despite interrupted flow upstream in the low gradient reach during 
the dry season. The Topanga watershed has a high water table relative to the creek and 
receives groundwater recharge in most reaches of the creek throughout the entire year (To- 
bias 2006). The Topanga Canyon rain gauge (#318) measured 11.73 cm of rainfall from 
October- December 2016, which was sufficient to breach the lagoon on 21 December 2016, 
but the breach was not sustained. Another 6.91 cm fell before the lagoon breached again 
on 11 January 2017. With an additional 42.11 cm of rain between January and February, 
Topanga Lagoon remained breached for a total of 132 d. It was the longest breach period 
for a larger surface area lagoon in this study and included an estimated 30 d of passable 
connection. The berm closed on 16 May 2017 and remained closed for the remainder of 
the study period. 

Topanga Lagoon is constrained on all sides by fill slopes, and PCH bridge abutments, 
which laterally constrain the lagoon’s physical features. This caused quicker and longer 
breaches than might have historically occurred, since the lagoon is not able to expand and 
spread out naturally to take up a larger surface area and volume.* Topanga Lagoon con- 
sistently maintained surface water throughout the drought period and required the second 
~ smallest amount of rain to breach in the 2017 water year. Tidewater gobies (Eucyclogob- 
ius newberryi) are supported in this lagoon, and while in-stream vegetation is limited to a 
small patch on the east side of the lagoon downstream of PCH, these are two potential 
food sources for smolts. 

Restoration that removes constraints on the sides of the lagoon above and below PCH 
would allow more area for important wetland and transitional upland bank vegetation in 
addition to potential for sea level rise accommodation. This must be balanced by preserv- 
ing enough thalweg constraint to continue maintaining a breach. Careful analysis of breach 
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flows needed to connect to the ocean, while expanding the duration of suitable flows that 
support fish migration will be needed. During the 2017 water year, two anadromous adults 
entered the system and one of them is confirmed to have spawned (RCDSMM unpub- 
lished data). Recent genetic analysis indicates that anadromous genes in this population 
are currently present, retaining the potential to produce anadromous individuals despite 
years of limited passage opportunities and purely residential spawning (RCDSMM unpub- 
lished data). Restoration planning was initiated in 2004 with the preparation of required 
documents for replacing the PCH bridge, which continues to date. 

Removal of three fish passage barriers in Arroyo Sequit in 2016 restored access to four 
kilometers of high-quality spawning and rearing habitat for O. mykiss, which was utilized 
by two anadromous adults in January 2017. Arroyo Sequit lagoon saw variable condi- 
tions throughout the drought period between 2013 and 2016, with consistent pooling un- 
der the PCH bridge associated with a high water table in that area. This relatively small 
watershed (31 km*) receives no imported water and is sustained only by groundwater flows 
from the upper watershed. The lagoon pool was not connected either upstream to the dry 
creek channel or downstream to the ocean except for one day on 29 October 2015 fol- 
lowing a storm. The lagoon had ponded surface water during the winter of 2016, as well 
as in September-December 2016. After a combined 18.92 cm of rain fell from October- 
December 2016, the first breach occurred on 20 December 2016. The breach closed for 
a brief period from 15-17 January 2017 but re-connected following a storm and remained 
breached through a combined 40.51 cm of rainfall over the following three months. Arroyo 
Sequit’s breach closed on 9 April 2017, but the lagoon remained subject to input daily with 
high tide overwash through the end of June 2017. The lagoon was breached for a total of 
97 d with an estimated 51 d of passibility. Schools of stripped mullet (Mugil cephalus) and 
Topsmelt (A therinops affinis) were documented in the lagoon but were not able to survive 
as the depth decreased to less than 30 cm (RCDSMM unpublished data). There are no 
proposals to make any changes to the lagoon at this time. 

In Trancas, the grade and concrete pylons/walls at the PCH bridge severely constrain 
the lagoon, and a channelized portion of the creek just upstream of the lagoon likely con- 
tributes to its current tendency to funnel water to the west side of the lagoon, aiding in 
a quicker breach. The current configuration of the PCH bridge severely limits the abil- 
ity of the lagoon to move as the bridge pylons redirect flow in a way that has built up 
sediment on the eastern side of the lagoon. Surface area of Trancas lagoon was in the 
mid to high range in our study, but also sustained a long breach period. The Trancas 
Lagoon Restoration Feasibility Study (Dagit et al. 2015) analyzed lagoon restoration al- 
ternatives to better understand these relationships and found that there were significant 
trade-offs between increasing lagoon and wetted area versus increasing or maintaining 
breach frequency. A larger lagoon footprint decreased the potential for breach, but in 
the face of sea level rise, adequate lagoon size will be necessary to buffer shoreline re- 
cession. While there are currently no O. mykiss present in this system, there is historic 
documentation of O. mykiss caught in the 1980’s (CalTrout 2006), the upstream habitat is 
good to excellent and most of the watershed is owned and managed by the National Park 
Service. 

Trancas lagoon experienced brief breaches during the 2014-2015 (5 d) and 2015-2016 
(1 d) water years. From October-December 2016, 17.63 cm of rain fell before Trancas la- 
goon breached on 24 December 2016. The breach was sustained as other storms dropped 
another 14.86 cm of rain from 30 December 2016 through 12 January 2017. The lagoon 
closed on 12 January 2017 but reopened again 8 d later with a new storm on 19 January 
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2017. The breach then remained connected until 18 April 2017, at which time the lagoon 
closed for the remainder of the water year. Trancas lagoon was breached for a total of 110d 
and passable for an estimated 61 d. The berm continued to build vertically as ponded wa- 
ter remained within the lagoon and under the PCH bridge through the end of the water 
year. Grooming of the berm by LA County Department of Beaches and Harbor occurred 
immediately after the breach ended. 

Restoration of Trancas lagoon is in the planning stages associated with the pend- 
ing replacement and expansion of the PCH bridge. This process is complicated by both 
landowner reluctance on the north side of PCH and the proposed re-nourishment of Broad 
Beach along the ocean. This project could potentially expand the width of the beach berm 
to almost 60 m, approximately twice its current width. Finally, the upstream concrete trape- 
zoidal flood control channels identified as the keystone passage barrier to the approxi- 
mately 10 km of upstream habitat are being studied to potentially restore a natural bottom 
that could facilitate fish passage opportunities. 

At Big Sycamore lagoon, both the PCH bridge at the ocean and a pedestrian bridge 
over the upper end of the lagoon connected by fill slopes on both sides entirely constrain 
the lagoon area. Waves actively deposit sand under the PCH bridge, creating a vertical 
berm that inhibits connectivity. The water is thus forced upstream above the pedestrian 
bridge and into the mainstem of the creek. Substantial vegetation lines the banks above 
the pedestrian bridge but no vegetation is present either in the lagoon or on the banks 
of the lagoon between PCH and the pedestrian bridge. The extent of the berm is exacer- 
bated by the concrete walls flanking the PCH crossing. Waves are able to deposit sand, 
but removal processes are hindered by the southern wall. Big Sycamore was ranked sec- 
ond in total habitat score for the Santa Monica Mountains (CalTrout 2006), but since then 
the Springs Fire in 2013 devastated much of the upper watershed and subsequent mud- 
slides and sediment deposition further degraded the habitat (NPS pers communication). 
Regrowth and improvement has occurred, and the 2017 water year showed promise in Big 
Sycamore creek, but access upstream 19 km to the limit of anadromy was minimal. The 
lagoon filled and a large reach of the creek upstream remained wetted well past the rainy 
season. There are currently no O. mykiss in this watershed. 

Big Sycamore lagoon completely dried down during the 2013-2016 drought period, 
with occasional puddles appearing in the lagoon due to sporadic overwash or a small 
rain event. In December 2016, the Lechuza rain gauge recorded 13.13 cm of rainfall, but 
Big Sycamore lagoon remained disconnected from the ocean. It was not until another 
28.37 cm of rain fell over two storm events that the lagoon finally breached on 23 Jan- 
uary 2017. After 5 d of connection, the berm closed, but was reopened one week later with 
another series of rain events that totaled 24.38 cm over a 31-d span. Despite this large 
amount of rainfall in the watershed, Big Sycamore lagoon closed on 8 March 2017. Small 
rain events in March and April 2017 were unable to cause another breach; the berm re- 
mained intact for the remainder of the year. An estimated 39 d of breach were recorded, 
the minimum for our entire dataset, with only 30 d of passable connection, most often 
at times of high tide. Once closed on 8 March 2017, the berm began to substantially 
build up. 

The PCH-pedestrian bridge complex has been identified as the keystone passage barrier 
in this system and has recently been placed on the priority restoration list by Caltrans 
District 7 (Caltrans pers. communication). Efforts are underway to obtain funding for the 
feasibility study needed to examine overall watershed processes post fire and develop a 
restoration plan. 
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The restoration challenges at Zuma lagoon are exacerbated by over 20 fish passage bar- 
riers extending from the Zuma Beach driveway and Arizona crossing, upstream through 
numerous private properties to the public area owned by National Park Service.!? Plan- 
ning has been in progress since 2005 to replace the Arizona crossing with a bridge, and 
to restore passage under PCH, but implementation funding is still missing. NPS restored 
over five acres of tules and transitional upland vegetation on the east side of the lagoon in 
2010, but management for mosquito abatement resulted in creation of paths through the 
tules that has fragmented that habitat and attracted transients. 

Habitat quality in Zuma lagoon was documented prior to the drought. Despite 
the added natural wetland vegetation, high water temperatures, dissolved oxygen level 
limitations, and extensive eutrophication would make this a challenging habitat for 
O. mykiss (McLaughlin et al. 2014). Zuma lagoon almost completely dried down during 
the 2013-2016 drought period and did not breach at all. Only a small fraction of the lagoon 
area had surface water in October and November 2016. It took a combined 42.09 cm of 
rainfall from October-January for Zuma lagoon to breach on 20 January 2017. Two more 
large rain events, with a combined total of 29.49 cm of rainfall recorded, sustained the 
breach for a total of 41 d with 39 d of passibility into the lagoon, dependent on high tide 
conditions. Once rain events ceased, the lagoon closed on | March 2017 at very high-water 
capacity, and the 200 m wide berm seen at this location built back up rapidly. Grooming 
and beach manipulation by LA County Department of Beaches and Harbors, who main- 
tain the berm closure to provide lifeguard vehicle access across the beach, assisted in the 
rapid berm buildup. 

Almost all the upstream fish passage barriers in Solstice have been removed during the 
past decade by the city of Malibu and National Park Service. The final constraint is a Cal- 
trans box culvert under PCH that has been the subject of much restoration analysis for 
years; it creates a steep >1 m drop-off with no jump scour pool that precludes fish passage 
during much of the breach period. Final planning for a bridge to replace the culvert is in 
progress but there is no current timeline for implementation. While there is historic docu- 
mentation of O. mykiss in Solstice (Dagit et al. 2005; NMFS 2012), lagoon development 
and presence is currently ephemeral and constrained by PCH and houses. Sometimes there 
is kelp washed up to the pool, but otherwise there is no vegetation and habitat is marginal 
at best. Flow through the culvert rarely exceeds 5 cm depth, although riparian cover and 
suitable habitat is available from the upstream end of the culvert to the limit of anadromy 
approximately three kilometers upstream. 

Solstice lagoon experienced breaches in the 2014-2015 (7 d), and 2015-2016 (11 d) water 
years. During 2016-2017, the first breach occurred on 23 December 2016 after a com- 
bined 14.73 cm of rain fell from 1 October 2016 through the day of first breach. With an 
additional 49.02 cm of rain recorded over the 2017 rainy season, the breach remained con- 
nected until closure on 4 July 2017. Solstice lagoon remained open for the longest period 
out of all the locations in our study area but was completely impassable. When the breach 
was first initiated, there was a short window before enough sand had washed away that 
would have allowed passage from the ocean to the culvert, but not further upstream. Sub- 


am Kelley, E., E. Wallace, and M. Stoecker. 2013. Southern California steelhead recovery and barrier 
removal recommendations for Zuma Creek, Santa Monica Mountains National Recreation Area. Pre- 
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sequent sheet flow across the beach face and the smooth concrete bottom of the culvert 
also restricts passibility. Additional human interference was a constant occurrence at this 
location. Sandbags were found within the culvert and obvious trenches appeared dug into 
the sand to divert flow. 

Although Caltrout (2006) and NMFS (2012) recovery plan both identify Las Flores as 
having potential O. mykiss habitat, there is no historic record of fish in this system and 
the anthropogenic constraints at the beach pose significant barriers. Las Flores lagoon has 
an extremely constricted channel bordered by concrete walls on either side both up and 
downstream of the PCH bridge, which funnels flows to the beach. Due to the fire hose 
funneling effect of these constraints, it breached every year despite drought conditions in 
the 2013-2016 water years. These breaches were variable in length, lasting 4, 10, and 7 d, 
respectively. A combined 16.51 cm of rain fell from 1 October 2016 until the first breach 
occurred on 24 December 2016. The breach was sustained throughout the rest of the rainy 
season, with an additional 39.12 cm of rain recorded, until its closure on 22 April 2017. 
Total breach time in 2017 equaled 119 consecutive days with an estimated 44 d of passibility 
into the lagoon but the lack of connected surface flow in the upstream channel precluded 
any further movement. Berm buildup occurred over the month of May 2017, and at the 
end of the 2017 water year what remained of the lagoon sat dry, 1-3 m below the height of 
the completely closed off berm. 

In addition to the site-specific constraints found at each lagoon system, the global ef- 
fects of sea level rise also need to play a role in developing restoration scenarios. According 
to a Los Angeles County Public Beaches Sea Level Rise Vulnerability Assessment (Noble 
2009), most beaches in the western portion of the county are expected to be eroded under 
a high sea level rise scenario by 2040. Topanga beach is expected to erode completely and 
Zuma beach is projected to erode by ~70%. By 2100, Zuma beach and Topanga beach 
are expected to recede by 184 feet and 190 feet respectively. While the other lagoons in 
our study were not included in the LA County Assessment, we can infer that erosion 
and recession could be similar under high sea level rise conditions. Explicitly evaluating 
the ecological trade-offs between ensuring adequate connectivity and increasing transition 
zone habitat in the face of sea level rise, as well as examining the best way to develop a 
realistic, self-sustaining mosaic of ecosystem habitat types and services, is necessary but 
difficult (Southern California Wetlands Recovery Project 2018). Continued monitoring of 
lagoon volume over time can help inform management decisions by providing accurate 
knowledge of how a given lagoon responds to rainfall, and how it provides habitat to emi- 
grating smolts and anadromous adults returning to riparian headwaters. It will be critical 
for managers and restoration planners at all lagoon sites to balance the need for breach 
connectivity with lagoon footprint as restoration plans move forward. 


Conclusions 


The recovery of anadromous adult steelhead trout to headwaters of southern California 
creeks is contingent on open and passable lagoon conditions. Our analysis of lagoons in 
the Santa Monica Bay indicates that lagoon surface area correlated positively with the 
amount of rain needed for the first breach to occur but negatively with the number of 
days a breach persisted. Larger lagoons in our study area tended to open later and close 
more quickly than smaller lagoons. At the same time, larger lagoon footprints could offer 
sea level rise protection to estuarian vegetation and other wildlife that is critical to the 
overall health of the ecosystem (Southern California Wetlands Recovery Project 2018). As 
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observed during the study period, lagoons with smaller footprints have little potential to 
provide smolt rearing habitat and provide limited anadromous adult passage. At this point 
in time, all of the lagoons in our study region would be considered unsuitable for rearing 
smolts. While it has been continually observed that O. mykiss smolts do not linger for very 
long in estuaries!! (Welch et al. 2004), insufficient lagoon volume can also lead to increased 
predation risk by birds (Kelley 2008; Osterback et al. 2015). Lagoons with larger footprints 
have better potential to be utilized as smolt habitat with decreased predation risk since they 
retained volume the longest through the summer and into the next water year. Adequate 
depth played a role here, as larger surface area locations were deeper overall, providing 
habitat that is more protected than shallower areas. The problem is that these lagoons 
also tend to breach less frequently, restricting anadromous adults from entering to spawn, 
which can result in fewer fish overall (NMFS 2012). 

The last important part of restoring connectivity for O. mykiss is for fish to have access 
from the lagoons upstream to headwater spawning and rearing habitat. Dry down of lower 
reaches in these creeks poses a significant barrier for both adult and juvenile migration to 
and from the lagoon. As extensively documented in Topanga Creek, the stream channel 
dries down from 200 to ~1700 m every summer (Dagit et al. 2017), effectively rendering 
the lagoon habitat useless to any anadromous adults or smolts that failed to enter the creek 
or ocean during a short period of passable connection. This pattern was also observed in 
all the other watersheds. Water temperatures in the study area lagoons consistently reach 
stressful levels during the summer months and these already high temperatures increased 
over time (RCDSMM unpublished data). Additionally, instream vegetation, which pro- 
vides habitat for benthic macroinvertebrates and terrestrial insects, and provides a major 
food source for juveniles (Krug et al. 2012), is virtually nonexistent in every lagoon of our 
study area (RCDSMM unpublished data). 

The measurements recorded and analyzed in this study provide insight into the volu- 
metric characteristics that small creek archetypal lagoons display in southern California’s 
Mediterranean climate. With the Santa Monica Bay predicted to face more extreme condi- 
tions due to climate change and sea level rise, maintaining sufficient volume in the lagoons 
in our study area as well as managing them to support connectivity will become increas- 
ingly more important to the recovery of southern steelhead trout. Unfortunately, due to 
the small data set with only three points in time measured for each lagoon, robust statis- 
tical analysis was not possible to further illuminate the complex dynamics that need to be 
understood in order to develop self-sustaining long-term restoration plans. 

Passage limitations due to drought are a key limiting factor affecting adult steelhead 
in the Santa Monica Bay, and despite the ability to spawn multiple times anadromous 
adult movement is hindered by flashy peak flows characteristic of southern coastal streams 
(NMFS 2012). In both Topanga and Arroyo Sequit, anadromous adults were able to mi- 
grate upstream in 2017, despite the short window of opportunity for passage. They were 
not able to return to the ocean due to low flow conditions in the creeks. If climate pre- 
dictions for low flow conditions continue and further restrict migration opportunities, 
O. mykiss may have no choice but to remain in freshwater refugia pools. Resident O. 
mykiss are typically smaller and less fecund than ocean-going steelhead (NMFS 2012). The 
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implications of reduced connectivity and migration opportunities could decrease resiliency 
of southern O. mykiss by increasing the risk of local extirpations and restricting genetic 
flow among neighboring systems (NMFS 2016). As sea level rise is projected to inundate 
wetland habitat across southern California in the next 100 years (Thorne et al. 2018), now 
is an opportune time to set up management plans to ensure adequate estuarine and transi- 
tional habitat remains resilient in the face of climate change and anthropogenic stressors. 
This study provides insight into how the lagoons in the study area rebounded from drought 
conditions and subsequently evolved over a wet period followed by the hot and dry season. 

With the current population of southern steelhead trout at record lows, the importance 
of passage opportunities in the Santa Monica Bay are more important than ever. Increased 
fire risk and uncertain precipitation levels that have come along with climate change will 
be an important consideration for management in the future as populations have the po- 
tential to be completely extirpated during a single natural hazard. Enhancing the ability of 
southern steelhead trout to move not only within their own watershed, but also between 
watersheds, could be vital to preserving this population. Finding the balance between la- 
goon volume and breach potential will be critical. 
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Obituary 


Image: Claremont McKenna College 


Long-time SCAS Fellow and Professor of Biology Emeritus at Claremont McKenna 
College Daniel Guthrie died on July 1, 2019. 

Guthrie, 80, specialized in comparative anatomy, animal behavior, introductory biology 
and environmental science during his 48 years at the Claremont Colleges, and had a lasting 
impact on students interested in careers in biology, ecology and medicine. Guthrie was an 
active member of SCAS in which he served as secretary, treasurer and Bulletin editor, and 
whose contributions to the academy’s mission are unprecedented. 

Dan authored more than 70 papers and studies in biology, and led innumerable field trips 
both locally and abroad. Every summer the Guthrie family drove across the United States, 
camping, looking at birds and collecting fossils. One of his favorite yearly expeditions was 
to San Miguel Island, which is one of the most remote destinations amongst California’s 
Channel Islands. He also had an extensive collection of seaweed that he gathered from 
Alaska’s Aleutian Islands while birding on Attu Island. His love for learning never faded, 
which inspired his travels to more than 40 countries, including Madagascar, Mongolia and 
Belarus. 

Born in Terre Haute, Ind. and raised in Garden City, N.Y., Guthrie attended Amherst 
College, where he studied biology and played lacrosse. He earned a Master’s degree in bi- 
ology from Harvard and a doctorate from the University of Massachusetts. As a biologist, 
he salvaged and collected U.S. bird skeletons for use by museums and teachers. He donated 
his time and expertise to finding and identifying fossil bird bones from roads, excavations, 
and collections for local educational institutions. 

Even as his mobility worsened, Dan still got out of the house to watch birds and enter 
data into Cornell’s eBird database. His extensive collection of bird watching notebooks 
continued through 2019, which included observations of more than 7,000 bird species. 
Guthrie also loved soccer, attending games at Claremont College, and watching the FIFA 
Women’s World Cup right up to his passing. 

SCAS offers the sincerest condolences to Dan’s family: daughters Kate Poaster, 
Ruth Guthrie and Winnie Larson, and their spouses, Bob, Dave Overoye and Tom, 
respectively; grandchildren Lee, Ben, Bobby, Acacia, Sage, Jack, Teddy and Alex; and 
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great-grandchildren Aydan, Gavyn and Cypress; and to all whose lives and careers were 
touched by Dan’s generosity and influence. 


Dave Ginsburg, Ph.D. 


SCAS President 
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